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Due to increasing importance of glass fibre 
reinforced plastics in various industries, it is necessary 
to develop accurate and efficient stress analysis techniques 
for them. With this in mind, the present investigations 
have been carried out to improve the brittle coating method 
and the technique of transmission photo elasticity to a state 
where they can be easily applied to the problems of stress 
analysis involving orthotropic materials. 

The studies with brittle coating have been per- 
formed using an indigenous coating developed in the 
laboratory which permits an easy variation in the proper- 
ties required for its application on the polymer matrix 
surfaces of composite materials. The coating behaviour 
has been studied on uni directionally reinforced glass 
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fibre epoxy composite subjected to uniaxial and biaxial 
stress fields. It has been established that the cracks 
represent the directions of principal strains and not the 
directions of principal stresses and that the threshold 
strain for the coating remains unaffected by the ortho- 
tropy of the specimen material. . Details of a data analysis 
procedure have been worked out so as to readily obtain 
specimen stresses from crack patterns. 

Optically superior photoelastic orthotropic model 
materials have been developed with glass fibres and polyes- 
ter resin. A systematic method of assessing the suitability 
of heterogeneous materials to photoelastic applications has 
been developed through considerations of laws of optics. 

The newly developed model materials exhibit better trans- 
parency, higher resolution of closely spaced lines, negli- 
gible distortion of a directed beam of light and sharper 
isochromatic fringes compared to the materials reported in 
the literature. 

An assessment of existing theories of photo- 
orthotropic-elasticity shows that some of the theories are 
capable of predicting photoelastic response of orthotropic 
materials quite accurately. However, no theory permits 
evaluation of the individual strain or stress components 
or even the difference between their principal values from 
birefringence measured in a polariscope. The directions 
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of principal stresses or strains also can not be determined. 
An exact strain-optic lav; derived on the basis of theory of 
photoelasticity for anisotropic crystals also suffers from 
these limitations. 

In order to overcome the above limitations in the 
application of photoelasticity to the stress analysis of 
fibrous composites, an approximate strain-optic lav; has 
been derived by minimising error between approximate and 
exact laws. A method has been developed to estimate the 
error introduced by the approximations. Thus, it has been 
possible to establish, with a reasonable degree of accuracy, 
the difference between principal strains at every point in 
a model. The theoretical predictions regarding the error 
estimation have been verified through experiments conducted 
on orthotropic circular disc under diametral compression 
and on rectangular specimens under uniaxial tension. 

With the above developments, the technique of 
transmission photo elasticity can now be applied to solve 
problems of the stress analysis of orthotropic elasticity 
as easily as it is employed for isotropic materials. 



CHAPTER 1 


INTRODUCTION 

1.1 INTRODUCTION TO COMPOSITE MATERIALS 

Demands of materials imposed by todays advanced 
technologies have become so diverse and severe that they 
often can not be met by single component materials acting 
alone. It is frequently necessary to combine two or more 
materials into a composite to which each constituent not 
only contributes its share, but whose combined action 
transcends the sum of the individual properties and pro- 
vides new performance unattainable by the constituents 
acting alone. Many engineering applications impose require- 
ments that are best met and in many instances met only by 
composite materials. Different types of composite mate- 
rials are developed employing various techniques best 
suited for a particular application. The fibre reinforced 
materials are probably the most important class of compo- 
site materials. They are very widely used in structural 
applications particularly when the weight saving is a pre- 
mium because they offer high strength to weight and modulus 
to weight ratios, excellent corrosion resistance} forma- 
bility and mechanical properties which can be controlled 
to a significant degree. For an efficient utilisation of 
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these materials and their adoption to newer applications, 
it is necessary to develop accurate and efficient stress 
analysis techniques applicable to such materials. 

1.2 LITERATURE SURVEY AND SCOPE OF PRESENT WORK 

An important characteristic of the fibre reinfor- 
ced materials is that they have high strength and modulus 
in the direction of fibres and the properties in other 
directions are different from those in the fiber direction 
that is; the properties change with the direction. Macro- 
scopically fibrous composites can be idealised as ortho- 
tropic materials, a class of anisotropic materials having 
three mutually orthogonal axes of material symmetry. Ana- 
lytical solutions to the stress analysis problems involving 
anisotropic materials are extremely difficult. The general 
approach to the theory of anisotropic elasticity has been 
discussed by Lekhniskii Ell. Analytical solutions to some 
simplified problems with thin infinite bodies have been 
obtained by Lekhniskii [ 1 1 Savin [ 21 and others E 3, k - 1 • 
However, analytical solutions to more realistic problems 
are not available. An alternative approach to solve such 
problems appears to be numerical one. Out of the various 
numerical methods, finite element technique is considered 
generally a very useful and powerful tool in analysing 
structural components. Complex geometries and complicated 
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loadings can bo handled with ease. Solutions to classical 
and special problems are being developed [5, 6] using 
this technique. However, a greater confidence into the 
numerical solutions can be placed, only when they are well 
supported by experimental observations. In fact numerical 
methods and experimental techniques can be employed to 
supplement each other. They two together can provide more 
information about the problem being analysed than either 
method taken alone. An example of this has been provided, 
by Sanford and Beaubient 71 who analysed a complex part 
made of single phase isotropic material. 

There are a number of experimental methods 
available to solve the stress analysis problems involving 
isotropic materials. Some of the methods can be directly 
applied to anisotropic materials also whereas others 
require modifications. Electrical resistance strain gauges 
are widely used for measuring surface strains. They can 
be applied to isotropic or anisotropic materials with equal 
ease. They are efficient and accurate in measuring strains 
at a point. Recently the techniques in strain gauges have 
been further developed to measure internal strains. For 
example, the gauges can be embedded in the composite mate- 
rial to measure residual strains [ 8] . The strain gauges 
can be employed very efficiently to measure strains when 
the critical areas are known, however, their use becomes 



very cumbersome and expensive if complete strain distribu- 
tions over large areas are to be dete rmin ed. 

Moire technique is also based on the measurement 
of surface displacements and hence can be applied to aniso- 
tropic materials. This technique has been employed to 
obtain strain distributions around holes and other free 
edges in laminated composites [9 - 11 1. However, this 
technique requires more sophisticated instrumentations 
which are not available in many laboratories. Holographic 
technique has also been extended for stress analysis of 
anisotropic materials [ 12], But this technique requires 
even more complicated and expensive set-up so that it may 
not be adopted as a routine stress analysis method. 

The present investigations have been carried 
out to develop techniques of brittle coating and trans- 
mission photoelasticity to a state where they can be easily 
applied to the problems of stress analysis involving ortho- 
tropic materials. Brittle coating technique is a very 
useful method for locating critical areas in structures and 
machine components through preliminary stress analysis. 

The critical areas can then be analysed using one of the 
more accurate stress analysis methods. The technique 
of brittle coating for isotropic materials is fairly well 
developed. Cunningham and Yavorskyt 131 applied brittle 
coating to wood which is an orthotropic material. 

However, there are many problems associated with the 
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application of this technique to orthotropic materials and 
in the interpretation of experimental observations. Such 
problems have remained unexplored yet, probably because 
of the limited use of wood in engineering applications. 

Now with the development of fibrous composites, a fresh 
look at the brittle coating technique and associated prob- 
lems is necessary. The present work gives a lead in this 
direction. 

Photoelastic stress analyses of orthotropic 
materials are performed either by applying a birefringent 
coating on structures and machine components or by making 
models of transparent birefringent materials. Dally and 
Alfirevich t 14- 3 applied an isotropic birefringent coating 
on glass fibre epoxy composite to obtain stress distribu- 
tions around a hole. They encountered difficulties due to 
mismatch of Poisson 1 s ratio between the specimen material 
and the birefringent coating. The situation is further 
complicated because the Poisson’s ratio of the specimen 
material changes with direction. Pipes and Dalleytl^l 
carried out a theoretical analysis which could be used in 
the interpretation of birefringent-coa ting-isochroma tic 
data for laminated-composite-material applications. They, 
however, pointed out that the birefringent coating method 
is likely to produce greater inaccuracies near the free 
edges like holes, cutouts etc. and also in the areas of 
high strain-gradients. 
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In recent years several investigators have success- 
fully initiated work on extending the technique of trans- 
mission photoelasticity to the stress analysis of fibre 
reinforced materials. There are two primary requirements 
for such applications; one, the composite model material 
should possess appropriate optical properties such as 
transparency and birefringence, and two, an appropriate 
theory to translate the photoelastic measurements into 
desired stresses and strains. It is interesting to note 
that although first attempt to develop transparent compo- 
sites for photoelastic analysis seems to have been made by 
Horridge [16] in 1955? the real photoelastic investigations 
on fibrous composites started much later with the publica- 
tion of results of Pih and Knight [17] in 1969 . Since then 
many investigators developed composite model materials and 
used them in photoelastic investigations. However, no 
scientific approach was developed either to produce opti- 
cally superior model materials or to assess their suitabi- 
lity to sophisticated applications. This has been done in 
the present investigations. 

Several theories of photo-ortho tropic-elasticity 
have been proposed in the literature. Some of them have 
been recently reviewed and assessed by Agarwal and 
Chaturvedi 1 181 . Pih and Knight [ 17 ] were probably the 

first to propose a stress-optic law for orthotropic 
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materials, however, their approach was unrealistic that 
they proposed a single photoelastic constant for orthotro- 
pic materials. It was not pursued any further by other 
investigators. Sampson [191 developed the concept of Mohr 
circle of birefringence by simple analogy and introduced 
the concept of three photo-elastic constants for orthotro- 
pic materials which are to be determined experimentally. 

Once the three constants are known for a material, the 
birefringent response of the material can be predicted 
through the stress-optic law formulated by Sampson. He 
supported his theory by limited experimental results. 

Dally and Prabhaka ran [20], like Sampson [19 1 , 
assumed that the photoelastic response of orthotropic mate- 
rials can be predicted through three independent photoelas- 
tic constants. Their experimental results for uniaxial 
state of stress agreed with their theory and with that of 
Sampson. Dally and Prabhaka ran [20], using fundamental con- 
cepts of fibrous composite, developed a method of theoreti- 
cally predicting the three photoelastic constants from the 
known photoelastic and elastic constants of the constituents 
and the volume content of glass fibres. Accuracy of such 
predictions is greatly influenced by the accuracy with which 
stresses in the constituents can be evaluated. This 
becomes a limitation of the approach particularly because 
it is very difficult to predict transverse stresses 
accurately. Another limitation of their work had been 
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that they did not discuss the isoclinics or their 
interpretation with regard to orthotropic materials. 
Prabhakaran [21] in a later paper did discuss the interpre- 
tation of isoclinics in photo-orthotropic-elasticity. His 
experimental results agreed well with Sampson’ s theoretical 
predictions for isoclinics. 

Bert [ 22 1 modelled fibrous composites by an 
orthorhombic crystalline system and applied the theory of 
photoelasticity of crystals first proposed by Pockels [23] 
and later modified by Bhagavantam [2*4-1 . Bert showed that 
the concept of a Mohr circle of birefringence as intuitively 
proposed by Sampson is a direct result of tensorial nature 
of birefringence. More recently, Knight and Pih [2^1, 
again using Bhagavantam theory, derived general stress - 
optic law and then simplified it to a two dimensional case. 
Their formulation is general and accounts for initial bire- 
fringence also. However, their simplified equations are of 
the same form as those of Sampson and Bert. 

Pipes and Rose [261 proposed the use of strain- 
optic law in place of stress-optic law. They assumed that 
the isotropic strain-optic law (with single strain- 
fringe value)is applicable for low fibre volume fraction 
birefringent composites. They experimentally verified 
this assumption for material with low degree of orthotropy. 
Prabhakaran [271 proposed a strain-optic law analogous to 
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that they did not discuss the isoclinics or their 
interpretation with regard to orthotropic materials. 
Prabhakaran [21] in a later paper did discuss the interpre- 
tation of isoclinics in photo-orthotropic-elasticity. His 
experimental results agreed well with Sampson’s theoretical 
predictions for isoclinics. 

Bert [ 22 1 modelled fibrous composites by an 
orthorhombic crystalline system and applied the theory of 
photoelasticity of crystals first proposed by Pockels [23] 
and later modified by Bhagavantam [ 2k] . Bert showed that 
the concept of a Mohr circle of birefringence as intuitively 
proposed by Sampson is a direct result of tensorial nature 
of birefringence. More recently, Knight and Pih [251, 
again using Bhagavantam theory, derived general stress - 
optic law and then simplified it to a two dimensional case. 
Their formulation is general and accounts for initial bire- 
fringence also. However, their simplified equations are of 
the same form as those of Sampson and Bert. 

Pipes and Rose [261 proposed the use of strain- 
optic law in place of stress-optic law. They assumed that 
the isotropic strain-optic law (with single strain- 
fringe value )is applicable for low fibre volume fraction 
birefringent composites. They experimentally verified 
this assumption for material with low degree of orthotropy. 
Prabhakaran [271 proposed a strain-optic law analogous to 
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stress-optic laws discussed in the preceding paragraphs. 

In the same paper, Prabhakaran has proposed a simplified 
strain-optic law with the single strain-fringe value which 
is taken as the arithmatic average of the three principal- 
fringe values. The approaches of Pipes and Rose and of 
Prabhakaran are really some kind of ad hoc approaches having 
no fundamental basis of material behaviour. Moreover, no 
added advantages have been gained by the use of strain- 
optic laws over the stress-optic laws. 

Besides what has been discussed above, photo- 
elastic investigations on composite materials have been 
carried out in other directions also. For example, the 
effects of residual birefringence t 28 - 30l and the study of 
stress-waves propagation in fibrous composites [ 31, 32]have 
also been reported in the literature. 

From the preceding discussion on the theories of 
photo-orthotropic-elasticity, it is quite clear that the 
photoelastic response of orthotropic composites can be 
accurately predicted through three independent photoelastic 
constants. It is, however, not clear if the existing 
stress-optic and strain-optic laws are adequate to inter- 
pret the photoelastic observations and translate them 
into the desired stresses or strains. Experiments were 
conducted on the superior unidirectional and bidirectional 
materials to assess the suitability of existing theories 
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of photo-orthotropic-elasticity. The results have been 
described in Chapter 4. Limitations of the existing stress 
optic laws have been pointed out. An exact strain-optic 
law has been derived using well established theory of 
photoelasticity of anisotropic crystals. An approximate 
strain-optic law has been derived in Chapter 5 from the 
exact strain-optic law. The approximations have been made 
to ensure minimisation of the error. Finally a method has 
been developed to estimate the error with a reasonable 
degree of accuracy. Experiments have been performed to 
verify the hypothesis of error estimation. 


CHAPTER 2 


BRITTLE COATING STUDIES ON EIBROUS COMPOSITES 
2.1 INTRODUCTION 

Brittle coating method is a very useful stress 
analysis technique which can be applied on actual structu- 
ral components. This avoids the preparation of models and 
conversion of results from model to prototype. Brittle 
coatings are widely applied for locating critical areas in 
structures and machine components through preliminary 
stress analysis. Its technology and the method of inter- 
pretation of results are well understood for the case when 
it is applied on isotropic materials. However, not much 
work seems to have been done for its application on ortho- 
tropic materials which may present many more problems. 

Cunningham and Yavorsky [13 1 applied the brittle 
lacquer on wood which is an anisotropic material. They 
performed tests under uniaxial state of stress and indicated 
that the coating cracks obtained on wood- specimens form 
trajectories of principal strains and not principal stresses. 
A transformation is needed to obtain the stress pattern in 
the specimen material. They had envisaged a future explo- 
ration of associated problems such as the relations between 
the law of failure of isotropic brittle coating and the 
strains and stresses in the anisotropic specimen. However, 
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such problems have remained unexplored yet probably 
because of the limited use of wood in engineering applica- 
tions. Now with the development of fibrous composites, 
a fresh look at the brittle coating technique and associa- 
ted problems is necessary. 

The present investigations were aimed at 
studying the behaviour of brittle coating on orthotropic 
materials and also studying the technological and inter- 
pretational problems. The coating has been applied to the 
uni directionally reinforced glass-fibre epoxy composite 
and its behaviour has been studied under uniaxial and 
biaxial stress fields. Differences arising due to the 
specimens being made of an orthotropic material rather 
than an isotropic material have been discussed. A data 
analysis procedure has been worked out so as to readily 
obtain specimen stresses from crack pattern. 

2.2 EXPERIMENTAL PROCEDURE 

Since commercial brittle coatings are not 
available in India, it was decided to develop a coating 
in the laboratory on the lines suggested by Ram Chand 
I 33 3 . Moreover, the indigenously developed coating 
permits an easy variation in the properties which is 
very helpful, almost essential, when the coating has 
to be applied on different materials such as epoxy or 
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polyester matrices used in the composites. Following are 
the basic constituents of a brittle coating; 

(i) Base Material; It is the most important constituent 
of the coating and significantly affects its ultimate 
properties. Zinc resinate has been used as the base 
material. Naturally occuring rosin and zinc oxide 
have been used for its preparation. The rosin is 
crystalline and very brittle in nature. It has a 
specific gravity of about 1.1 and melting point 

84- °C. Its natural colour varies from pale to dark 
yellow, which, however, does not affect the colour 
of the resulting coating. 

(ii) Solvent and Thinner; For brittle coating applica- 
tions, besides being a good solvent for the base 
material, the solvent should evaporate also reason- 
ably fast when exposed to atmosphere. Carbon 

di- sulphide has been used as the solvent and thinner 
in the present coating. However, it must be men- 
tioned here that carbon di-sulphide is an extremely 
hazardous solvent. It can easily ignite and its 
vapour, if inhailed, can cause health hazard. 
Therefore, extreme precautions should be taken in 
its use as described in reference [343. 
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(iii) Plasticiser? The plasticiser is added to control 
the strain- sensitivity and brittleness of the 
coating. Unplasticised coating may tend to "craze" 
(random cracks formed due to residual stresses) 
whereas, on the other hand, excessive plasticiser 
will reduce the brittleness so much that cracks 
once formed will close upon the release of the strain 
which caused the crack formation. In general, the 
plasticiser adversely affects adhesion characteris- 
tics of the coating. Dibutylphthalate is a common 
plasticiser. In the initial stages of the present 
coating development, it was used. However, it was 
affecting the adhesion characteristics too adversely 
and hence had to be dispensed with. The present 
unplasticised coating did not exhibit any tendency 
to craze . 

2.2.1 Coating Preparation 

The coating is prepared by first cooking 100 gm 
rosin with 1.6 gm zinc oxide at a temperature of 220 + 5 ° G 
for about 30 minutes. Continuous stirring is necessary 
during cooking to achieve a thorough mixing of constituents. 
The resulting resinate, a dark brown semi-solid material, 
is allowed to cool down to room temperature and then 
dissolved in 100 cc. solvent, that is, carbon disulphide. 
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The solution is filtered to remove unreacted elements. 
Additional 100 cc of carbon di- sulphide is added to the 
filtrate. The resulting dark brown solution is the desired 
brittle coating ready to be sprayed on to the specimens. 

It may, however, be pointed out that before finally arriving 
at the above procedure, several other procedures were tried. 
Relative proportions of the rosin and zinc-oxide were chan- 
ged as was curing time and temperature. A plasticiser was 
also used in the initial stages of development, 

2.2.2 Specimen Material 

Brittle coating studies were performed on 
commonly used fibrous composites. Uni directionally reinfor- 
ced glass-fibre epoxy composite was selected as the specimen 
material. The specimen material which was also fabricated 
in the laboratory exhibited the following average elastic 
properties in a tension tests 

Longitudinal • modulus, = 10005.5 MN/m 2 

Transverse modulus , E^ = 3835.5 MIT/m 2 

Shear modulus , & LT = 1481.0 MN/m 2 

Major Poisson' s ratio , v ^rj, = 0.38 

For studying the behaviour of coating under uniaxial state 
of stress, rectangular cantilever beam specimens of size 
125 x 20 x 3.5 mm shown in Figure 2.1a, were machined. 
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The fibre orientation, 0, was kept at 0, 30, 4-5, 75 and 

90 degree’ in different specimens. The behaviour of coat- 
ing under biaxial stress-field was studied through circular 
disc specimen of 100 mm diameter. The disc was loaded in 
diametral compression as shown in Figure 2.1b. Angle bet- 
ween the loading axis and the fibres was kept at 1 5 , 4-5, and 
60 degree in different tests. 

2.2.3 Application of Coating to the Specimen 

The coating was sprayed on the epoxy composite 
specimens degreased by acetone and further cleaned by 
carbon di-sulphide. An air compressor, which maintained a 
gauge pressure of 0.14- MPa + 10 percent, was used for 
spraying. In order to obtain a uniform thickness of coat- 
ing the spray gun was slowly moved along the length of the 
specimen and the gun was kept at a height of about 10 cm 
from the specimen surface. The number of passes of the 
spray were adjusted to obtain a coating thickness of about 
0.4- mm. 

The coating was cured in a temperature controlled 
oven at 120 °C for 2 hours and then the coated specimens 
were allowed to cool in the oven itself so that a slow 
rate of cooling is maintained. Rapid cooling outside the 
oven may cause thermal erazing of the coating. The coat- 
ing thus obtained exhibits satisfactory properties for 




18 


carrying out brittle coating stress analysis. Properties 
in question include adequate sensitivity, good adhesion to 
the specimen surface and clearly visible cracks. The 
cracks produced during loading remain open after removing 
the load. However, if there is a significant rise in the 
room temperature, the cracks tend to close partially or 
completely. The tests were conducted under conditions of 
controlled temperature and humidity at 25 °C and 4-0 per- 
cent relative humidity. 

2.2.4 Calibration of Coating 

Calibration of coating consists of finding out 
the threshold strain for the coating. The threshold strain 
is defined as the minimum strain which will produce cracks 
in the coating. It was obtained by applying the coating 
to a specimen with 0 - degree fibre orientation. The 
specimen was loaded as a cantilever beam and a known 
deflection (4- mm) was given to its free end. Position of 
the crack nearest to the free end was noted. The strain 
at this point gives the threshold strain for the coating. 
Magnitude of the strain was obtained experimentally by 
strain gauge measurements. The electrical resistance 
strain gauges were mounted at the desired position on 
the same specimen after removing the coating. The same 
known deflection was given to the free end and the strain 
was measured by a strain indicator. The magnitude of the 




Left to right: Fibre orientations 5, 30,45 
75 and 90 degree. 






22 

specimens are orthotropic and the directions of principal 


stresses and strains are, in general, different. The 
direction of major principal stress is the x - direction 
(Figure 2.3). The directions of principal strains x-rere 
experimentally obtained by strain gauge measurements. The 
angle, 3 , between the directions of principal stress and 
principal strain in the specimen has been calculated 
theoretically also as follows. A uniaxial stress, cj x , 
produces the following normal and shear strains in the 
specimen; 


e = 
x E x 

v- 

e y ~ ” xy ^ 
a x 

\j = ~ m x 


( 2 . 1 ) 


where EL, v and m„ are the elastic constants in 

x 7 xy x 

the x - y directions and are related to E^, E^, G^ and 
through equations given in Appendix A. It may be noted that 
the shear coupling coefficient, m x , is, in general, not 
zero. Due to the presence of, m x , directions of principal 
strains differ from the directions of principal stresses in 
the specimen. The angle, B , for the stresses and strains 
given by Eq. (2.1) can be shown to be 
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( 2 . 2 ) 
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The angles, g , calculated from Eq. (2.2) for 
different fibre orientation, and experimentally obtained 
from strain gauge measurements, have been compared, in 
Table 2.1, with the angles, ct , measured with the help of 
an optical projector (Baush and Lomb) from the brittle 
coating crack patterns on the corresponding specimens 
(The crack patterns are shown in Figure 2.2 and the angles, 
a , and 3, defined in Figure 2.3). It can be noted 
from Table 2.1 that the angles, a , and 6, have quite com- 
parable values. It may, therefore, be concluded that the 
brittle coating cracks indicate the directions of principal 
strains in the specimen and not the directions of principal 
stresses. This conclusion is of particular significance 
for the orthotropic materials where the two directions are 
different. In case of isotropic specimen materials, this 
conclusion is redundant. 

It is interesting to examine the effect of ortho- 
tropy of specimen material through a partially reinforced 
specimen. Figure 2.4- shows crack pattern on a specimen 
which contains inclined unidirectional fibres in part of 
its length while the rest was purely resin. The reinfor- 
ced and unreinforced regions have been shown separated by 
a dotted line and the direction of fibres is indicated by 
a double arrow. The cantilever specimen was deflected by 
an end load. It is seen from Figure 2.4- that the cracks 
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Table 2.1 


Results for Coating Behaviour under 
Uniaxial Stress 
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m the unremforced region are perpendicular to maximum 
principal stress whereas in the reinforced region it is 
not so although uniaxial state of stress exists everywhere 
m the specimen . 

In each of the tears conducted on the specimens 
with different fibre orientation, the minimum strain which 
causes the coating cracks has also been obtained from 
strain gauge measurements as well as from theoretical cal- 
culations. The valuesj thus obtained are within 4- percent 
of the threshold strain obtained earlier during coating 
calibration. This indicates that the threshold strain or 
the coating sensitivity is not significantly influenced by 
the degree of mismatch between the Poisson's ratios of 
specimen material and the coating. This is also important 
because, m an orthotropic material, Poisson's ratio changes 
with the direction, 

2.3.2 Coating Behaviour Under Biaxial Stresses 

The cs rcular disc specimens, coated with brittle 
coating, were compressed by two equal and diametrically 
opposite forces as shown schematically in Figure 2.1b. 
Diametral compression produces a biaxial state of stress 
at an interior point m the disc. Photographs of result- 
ing crack patterns tor fibre orientations of 1 5, 45, anc ^ 

60 degree are shown m Figure 2,5* The crack pattern is 
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clear throughout the disc except m the close vicinity of 
the points of application of load# Due to very high compre- 
ssive stresses m this region, the coating peels-off the 
surface of the specimen# 

For an interpretation of the crack patterns, it 
is required that the directions of principal stresses and 
strains be known at every point m the disc. Okubo [ If] 
reported an analytical solution of the problem of ortho- 
tropic circular disc under diametral compression. While 
trying to compute the stresses m the present orthotropic 
disc, some obvious errors were noticed m his formulation, 
the most important being an inconsistency m coordinate 
transformation. Therefore, even after considerable efforts, 
it was not possible to use his results with the present 
findings. One may probably develop a consistent solution 
to this problem using Okubo’ s approach. However, this was 
not done by the author for obvious reasons. Photoelastic 
analysis of ortho tropic disc under diametral compression 
has been attempted [21] . Isoclimc and isochroma tic 
fringe patterns have been obtained but it has not been so 
far possible to establish the directions of principal 
stresses or strains. It was, therefore, decided to esta- 
blish the directions of principal strains at three diffe- 
rent points in the disc using strain gauge measurements. 

The three points considered are shown m Figure 2.6. They 
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were kept the same for all orientations of the disc. The 
strain gauge measurements provide not only the directions 
of principal strains but also their magnitude. Therefore, 
the normal and shear stresses along the directions of 
principal strains can be found using the following stress- 
strain relations [351 . 

cr o x 

e = -£ _ v -S „ m JES 

P E p qp E q m p ^ 



Where the directions p and q refer to the 
directions of principal strains so that y pq vanishes. 

The elastic properties 3 p , E qJ G pq , v pq , v^, m_ p and m q 
of the specimen m the p - q directions can be obtained 
using transformation relations given m Appendix A. 

Eq. (2.3) can be solved for 0 , o q and t from which 
magnitude and directions of principal stresses can be easily 
estab] ished* 

The results of above measurements and analysis 
have been presented m Table 2,2 by giving angles which 
the maximum principal strain and maximum principal stress 
at a point make with the load axis. Also given m Table 
2,2 are the angles which the tangents to the cracks at 
the corresponding points malic with the direction 
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Table 2*2 s Re sails for Coating Behaviour Under Biaxial 
Stresses 


f 

Fibre orientation \ 

angle, ©, degree | 

! 

. L 

Point under 
consideration 

■ - . . 

a * 

CT 

degree 
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t rv ** 

i °e 
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degree 
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[ Jc 

} a crack 
i degree 


1 

85.7 

77.6 

76.5 

15 

2 

88.8 

85.7 

85.0 


3 

76.7 

66,2 

66.0 


1 

81.6 

89.0 

89-5 

45 

2 

78.0 

88.4 

88.7 


3 

7 0.7 

80.5 

83.5 


1 

88.5 

82.0 

81.5 

60 

2 

86.4 

83.6 

83.0 


3 

77.5 

82.5 

86.0 


= Angle between maximum principal stress and 
the load axis. 


** ct £ = Angle between maximum principal strain and 
the load axis. 


a = Angle between tangent to the crack and 
crack direction perpendicular to the load axis. 
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perpendicular to the loading axis. It can be clearly seen 
from these results that the cracks represent the directions 
of principal strains rather than the directions of princip®^ 
stresses. This observation is very significant for future 
application of brittle coating to the problems of stress 
analysis involving orthotropic materials. 

2 . b ANAL rat S TECHNI QUE 

It has been established m the preceding section 
that the coating cracks, produced under both uniaxial and. 
biaxial states of stresses, represent the common directions 
of principal strains and principal! stresses in the isotro- 
pic coating and represent the coincident principal strain 
directions in the orthotropic material underneath. Analysis 
procedure, based on these facts and aimed at evaluating 
specimen stresses m terms of coating stresses, has been 
developed m the following paragraphs. 

Lot p, q be directions, normal and tangential to 
crack at a point so that they represent directions of 
principal strains m coating as well as m the specimen 
and. also directions of principal stresses m coating. r ^ Le 
coating strains are related to costing stresses by isotro- 
pic stress-strain relations as follows s 


33 


c 


e 

P 



c 


o 

p 


C 0 

v cr 


q 



( 2 » 


Specimen stresses and strains are governed by orthotropic 
stress- strain relations - 


P 


Y 


pq 


JL 


E" 


P 


a 


q 


v 

pq 


qp 


CT S 

P 
,s 


e: 


q 


+ 


E 


P 


e; 


o = - 


P 


c s 

_J 2 _ 

Ti 1 ® 


m„ 


_PQ 


4 


- m 


s 

T 

J3SL 

4 


ni 


q 


a 3 

m ^ 1 .11 

s 


E. 


L 


( 2 . 5 ) 


G s 

pq 


Superscripts c and s denote coating and speci- 
men respectively. Since coabmg is usually very thin relative 
to Hie thickness of specimenj at is valid to assume that 
surface strains occurmg on the specimen are transmitted to 
the coating without magnifi cation or attenuation so that 

c s 

e = e 

P P 

( 2 . 6 ) 

c s 

e = e 

q q 

Tt is also reasonable to assume that a state of plane 
stress exists m both the coating and the boundary layer 
of the specimen. 
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Substitution of Eqs . (2.4) and (2.5) into Eq. 
(2.6) and solving for specimen stresses yields? 
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(A - B v c ) a c p + (B - A v C ) a G q 
AC - B 2 

C / -p C \ c 

a a+ (B~Cy ) q p 


(C - B v°) 


AC - B‘ 


(2.7) 
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It will be possible to calculate the specimen stresses from 

c c 

above equations provided the coating stresses and ° ^ 

are known at a point for a given load. Hie coating stre- 
sses can be accurately estimated at a point only when 
the crack appears there. It is assumed that at the minimum 
load aL which a crack appears, the coating stresses 


a 


and a q at the point of appearance of crack, are 
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related to the ultimate strengths of the coating through a 
failure law. For example, when the state of stress in the 
specimen is unidirectional, a crack arrears when the 
stress normal to the crack exceeds the ultimate tensile 
strength of the coating. Under the state of biaxial 
stresses, coating stresses can be estimated for any given 
load only after the cracks at the point under considera- 
tion have been observed m the two perpendicular directions. 
The estimation of coating stresses will require incorpora- 
tion of a suitable failure theory and the knowledge of 
ultimate strengths of the coating. Some of the failure 
theories have been discussed m Appendix B, Also given m 
the Appendix B is the method, of evaluating coating ultimate 
strength. The method of evaluating and under 

biaxial state of stress can be better illustrated through 
the following example m which Mohr theory as discussed in 
Appendix B has been employed to predict failure of the 
coating. 


c 

Let a be the tensile stress at the minimum 
P 

load P at which a crack first appeared at a point. To 

0 

obtain information about the cr^ the following two cases 

should be considered: 

c c 

Case 1 ( a > 0) : When cr is tensile, crack in 

— ~3 1 c 

the direction perpendicular to will appear for an 
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increased load P/ctP (c? < 1 ). According to Mohr theory 
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or a , 
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Hie specimen stresses 
through Eq. (2.7) by 
coating stresses, 
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Case 2: ° < 
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at load F 


( 2 . 8 ) 


at ] oad P can now he obtained 
substituting the above values of 

0 

s When a is compressive, the 


cracks in the perpendicular direction will not appear by 

increasing the load beyond P. In this case the magnitude 
c 

of is obtained by the technique known as relaxation 

technique t 36 , 37 1 . In the relaxation technique, an 
uncoated specimen is subjected to a known load. The coat- 
ing is applied to the specimen under load. Thus during 

c 

unloading of tne specimen coating stress a ^ will be 

c 

tensile and a will be compressive and the magnitude of 

ir 

stresses will be proportional to the amount of unloading. 

Under these condition 3 the cracks will appear m tho 

c 

direction perpendicular to a when the load has been 
P 

reduced by Y P, An estimation of stresses con be made 


as follows. 


37 


Let the coating stresses at load P during 
observations of first cracks be s 
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= a 
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-g P a 
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(2.9) 


Therefore the coating stresses in the relaxation procedure 
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when the load has been reduced by YP will be 
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where B and Y are both positive. The stress conditions 
indicated by Eqs. (2.9) and (2.10) will satisfy the failure 
conditions (B.2) of the Mohr theory oO that 
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Tlius the coating stresses a l load P ares 
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CHAPTER 3 


DEVELOPMENT AND CHARACTERISATION OF OPTICALLY SUPERIOR 
PHOTO ELASTIC ORTHO TROPIC MODEL MATERIALS 

3.1 INTRODUCTION 

The first reauireraent of applying the transmission 
photoelastic technicrue to the stress analysis of orthotropic 
materials is to have a suitable model material. A photo- 
elastic model material should be transparent, optically 
homogeneous and biref nngent . Glass, epoxies and polyester, 
which are constituents of many common orthotropic composites, 
are all photoelastic materials although they differ from 
one another in their optical properties. Therefore, glass 
fibre reinforced epoxy or polyester composites are natural 
choice for photoelastic orthotropic model materials. 

First attempt to prepare photoelastic comoosite 
material seems to have been made by Horndge f 16 1 who used 
a resin with refractive index close to that of glass fibres. 
Pih and Knight [ 17 ] produced orthotropic model materials 
with an enoxy resin and E - glass fibres by employing wet 
filament winding under vacuum, Sampson [ 19 1 used a dry 
filament binding of HTS glass roving followed by immersion, 
under vacuum, m an epoxy resin bath. Dally and 
Prabhakaran f 20 1 produced large sheets of E~ glass fibre/ 
polyester composites with unidirectional and bidirectional 
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reinforcement. The matrix material was blended with 30 
percent styrene monomer to improve its refractive index 
to match it with that of fibres. The model materials used 
by different investigators may be characterised as { suffi- 
ciently transparent’ , They have been used to obtain 
isoclinic and isochromatic fringes. However, the fringes 
are not sharp enough to be used for accurate stress analysis 
purposes. Some of them have not been found adequate m 
sophisticated applications such as dynamic photoelasticity 
[31 ], An improvement in the optical properties of model 
materials is definitely needed. This has been achieved in 
the present investigations first by understanding the effects 
of irihomogeneity through the laws of optics and then by 
developing scientific methods of assessing suitability of 
model materials to the photoelastic applications, 

3.2 OPTICAL CONSIDERATIONS FOR A MODEL MATERIAL 

A photoelastic model material should be transparent, 
optically homogeneous and birefnngent . Glass, epoxies and 
polyesters are all photoelastic model materials. However, 
when the glass fibres arc embedded m an epoxy or polyester, 
the resulting material is not as transparent as its consti- 
tuents. Decrease m transparency occurs due to scattering 
of light at the fibre matrix interfaces. Transparency is 
very significantly affected also by entrapped air m the 
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another effect which 16 more detrimental than the lost, of 
transparency Tor photo elastic applications, Refraction at 
the interface of such j uhomogonoi tics causes the path of 
light to deviate from its original straight line path. 
Moreover, in glass fibio reinforced composites the inter- 
faces are not straight so that non uniform deviation of 
light occurs from different points at the interface, Hie 
effect of refraction may be illustrated by considering the 
path of light, according to the laws of geometrical optics, 
through homogeneous and inhomogeneous materials as follows. 

Consider that a parallel beam of light is incident 
on a plane surface of a solid transparent medium (homoge- 
neous or heterogeneous) and that the surface, out of vhich 
the light emerges, is also plane and parallel to the surface 
of the beam incidence. In case of a homogeneous medium 
(Figure 3.1 a) the emergent beam of light will also be a 
parallel beam travelling m the direction parallel to the 
incident beam. The same is true of a medium which consists 
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of two different materials having different refractive 
indices cut uith the interface plane and parallel to the 
surfaces of light incidence and emergence (Figure 3*1 b) . 

In glass fibre reinforced materials, the glass phase is 
uniformly dispersed m the form of thin circular cylinders. 
The cross section of a composite containing a single fibre 
may he represented as shown m Figure 3*1 c. The parallel 
beam of light is shown to make a normal incidence on the 
material surface. It will continue to travel m the direc- 
tion of incident beam until it meets the fibre-matrix, 
interface. At the interface, due to change m refractive 
index, the incident ray deviates. The amount of deviation 
of a particular light ray will be determined by the orien- 
tation of interface and the degree of mismatch of refractive 
indices of the two media. Thus the parallel rays get 
deflected by unequal angles from different points of fibre- 
matrix interface resulting into diffusion of light beam. 
Conseauently, a finite-width incident light beam vill appear 
attenuated. The following simple experiment was performed 
to demonstrate the above effect. 

A model was prepared by embedding a parallel row 
of circular rods ±n a plate of an epoxy resin. The influe- 
nce of placing tins model in the path of a directed light 
beam has been examined. For this purpose, light caning 
out of a sJit source is passed through a set of two lenses 







so that the emerging beam is convergent one. A screen is 
placed at the focal plane oi the second len j so that a 
sharp image of the slit is formed on the screen. An 
arrangement of the optical bench is schematically shown m 
Figure 3 * 2 , When the model i d placer, between the second 
lens and the screen, a distortion in the image of the slit 
occurs. The direction of distortion changes with the orien 
tation of glass rods with respect to the slit. The images 
of slit have been shown m Figure 3*3 a to 3*3 c for 0, 4-5 
and 90 degree orientations and compared with the image 
without model (Figure 3*3 d) . It may be observed that the 
distortions m the image are always perpendicular to the 
axis of the glass rods. When the model is replaced by a 
plate of the epoxy resin without glass rods, no distoxtion 
of the image was observed. This simple experiment demons- 
trates that a glass fibre reinforced material may cause 
smearing of directed beam of light as schematically illus- 
trated m Figure 3*1 c. 

The affects of placing a glass fibre reinforced 
composite model m the path of light can now be summarised 
as follows i 

i) A parallel beam of light appears smeared m the 
direction perpendicular bo the fibres. 
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epoxy model. Angle of glass rods with the slit 
(a) 0°, (b) 45° and (c) 90°, and (d) image through 
unreinforced epoxy model 






li ) Sharp boundaries of an object, when viewed 

through the model, appear smeared, again, m the 
direction pernendicular to the fibres. As a 
consenuence very closely spaced straight parallel 
lines may not appear distinct. 

m) White light -passing through the model gets split 
m its colours and exhibits a double spectrum 
symmetrical about the centre line. This is caused 
because the refractive index depends upon the wave- 
length of colour of light. Thus, the rays of 
different colours deviate by unequal angles from 
the same point of the interface. 

The above effects are undesirable for photo- 
elastic applications. The first two cause a reduction m 
the sharpness of isoclinic and isochromatic fringes. The 
third one will interfere m l he analysis if white light is 
used. It is, therefore, auibe clear that the effect of 
mismatch of refractive indices is not only to reduce the 
transparency but to produce above undesirable effects also. 

It is, therefore, suggested that the matching of refractive 
indices be judged by the absence of above effects and not 
just by visual inspection for t ran soar ency. The first 
two effects have been utilised in the present investiga- 
tions to demonstrate optical superiority of the model 
materials. The third effect has been found to be very 
convenient for auick visual inspection of the material. 
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3.3 fabrication OF HATEFUL 

The uwo most important conditions for obtaining 
optically superior composite model materials ares 

i. the refractive indices of the fibres and the 
matrix material should, be matched as closely as 
possible and 

ll. the entrapped air in the composite should be 
minimised. 

The model materials developed in the present study 
are made of E - glass fibres (a product of Fibreglass 
Pilkmgton Ltd., India) and polyester resin. The polyester 
resin of a low viscosity and clear casting grade (a product 
of Glafpster Industries, Surat, India) has been used with 
0.4 per cent cobalt octate as accelerator and 1.0 percent 
methyl ethyl ketone peroxide (MEKP) as a catlyst. The 
cured resin exhibits a refractive index of 1.5532 m the 
white light. This refractive index can be increased by 
adding di vinyl benzene and decreased by adding di butyl ph- 
thalate to the resin. The refractive index of glass fibres 
has been quoted to be 1.549 + 0,003 by the manufacturer. 
However, the value has not been found to be very reliable 
one. Moreover, the refractive index may get affected by 
the coupling agents and coatings applied to the surface 
of the fibres. Consequently, refractive index of glass 
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fibres m a roving is different from that of glass cloth. 
After trying different percentages of above additives it 
was found that an addition of 1.9 per cent di butyl phthalate 
by weight of polyester resin accurately matches its refrac- 
tive index with that of the glass fibres from roving which 
were used for unidirectional composites. Additions of 5.0 
per cent dibutylphthalate and 2.0 per cent divinyl benzene 
are needed to match the refractive index of resin with that 
of the glass cloth which was used for bidirectional compo- 
site, It may be pointed out here that the refractive index 
of the liquid resin is different from that of the cured 
resin. Hence, matching of the refractive indices has to 
be achieved when the resin is completely cured. 

The composite plates were cast m a mould made 
of two plexig] ass plates separated by tough rubber spacers 
of uniform thickness. The glass fibres were would on a 
stiff, plane, rectangular frame of G.I. wire. The two 
ends of the fibre bundles were tied to the frame while 
maintaining tension m the fibres. The frame along with 
the fibres, gets embedded in the resin and is irretrievable. 
The test specimens have to be prepared from central por- 
tion of the cured plates away from the embedded frame. 

The resin and other chemicals were thoroughly mixed and 
poured into the mould after introducing the frame m it 
through the side opening. There are two primary sources 
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of air entrapped m the material prepared by this method. 
When the icsm and other chemicals are mixed and poured into 
the mould, entrapment of some air is unavoidable. The ali- 
en trapped during mixing i s removed by placing the mixture m 
a vacuum chamber for sufficient time prior to its pouring. 
After pouring, the mould is again placed m the vacuum 
chamber. Photograph in Figure shows a mould m the 
vacuum chamber. This removes air entrapped during pouring 
and also due to improper impregnation of fibres which is the 
second major source of entrapped air. Hie improper impreg- 
nation occurs because the resin is some- times unable to 
penetrate through the fibre bundles and wet all the fibres. 
To improve impregnation, a low viscosity resin was selected 
for the studies. Secondly, the fibre bundles were broken 
into smaller groups while winding them manually on the wire 
frame. This enhances resin penetration. Breaking of few 
fibres due to their handling m tins manner may slightly 
reduce the ultimate strength of the composite but ic does 
not affect the photoelastic investigations appreciably. 

The problem of incomplete wetting is more serious in glass 
cloth and twisted fibres. 

For making bidirectionally reinforced composite 
material, glass cloth was cut into pieces of proper size 
and were placed m a mould one by one. For better impreg- 
nation, a small quantity of resin was poured on each layer 
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Figure 3.4 A mould in the vacuum chamber to remove 

entrapped air. 





and brushed. Eight layers of cloth were used m a plate. 

The mould was placed in a vacuum chamber to remove entrapped 
air. 

In each case, the plates were cured m the mould 
at room temperature for 48 hours and post cured at 80 °C 
for 6 hours m a temperature controlled oven. For the pur- 
pose of studying optical properties of the materials, circular 
discs of 50 mm dia. were cut from each plate and machined. The 
discs were examined for optical properties discussed in the 
preceding section. The materials have been found to be 
adequate model materials. It may, however, be mentioned 
again that the proper chemical composition of the resin was 
arrived at after lot of trial and error. In the process, 
many compositions were tried anci the composite plates made 
but later discarded due to inadequate optical properties. 

Some of the better materials among the discarded ones have 
been compared with the model materials to illustrate latter’s 
optical superiority. 

3.4 OPTICAL SUITABILITY OF MODEL MATERIALS 

Optical properties of the unidirectional and 
bidirectional composite model materials have been examined 
to assess their suitability for photoelastic analysis and 
allied applications. The model materials have been tested 
for transparency, resolution of closely spaced lines, image 
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formation (or distortion) and photoelastic effect. The test 
procedures and results are discussed in the following para- 
graphs, 

3.^,1 Transparency- 

Transparency of model materials has been examined 
m the plane polarised monochromatic sodium light which is 
commonly employed for photoelastic investigations. The 
circular disc model was placed m a beam 01 convergent light 
as shown schematically m Figure 3.2. In each case, the 
setting of optical bench was kept the same and the model 
was placed at a fixed distance of 15 cm away from the screen. 
Intensities of the incident and emergent light were measured 
by a photometer (Model 520 M, Photovolt Corporation, New 
York) by placing its sensing element at the centre of screen 
without and with the model placed in position respectively. 
Thus, the transmission ratio defined as the ratio of inten- 
sity of emergent light to the intensity of incident light 
has been obtained in each case. 

The transmission ratios have been lound to be 
0.58 for the unidirectional model material and 0.55 for 
the bidirectional material. The transmission ratio for 
the superior materials is not influenced by the orienta- 
tion of fibres with respect to the slit of the light 
source, A slightly lower value of transmission ratio for 
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bi directional composite indicates that it is slightly less 
transparent because it has higher volume fraction of glass 
fibres. The transmission ratios of the two composite model 
materials are lower than the transmission ratio of the 
unreinforced polyester disc of equal thickness which was 
found to be 0,8. However, both the materials are suite 
adequately transparent. Another model of glass fibre polyes- 
ter composite with only marginally different resin composi- 
tion (consequently, slightly bad machmg of refractive 
indices) showed a lower transmission ratio of 0.50 when 
fibres were perpendicular to the slit and further reduced 
to 0.4-0 for fibres parallel to slit. A glass fibre epoxy 
composite, which looked adequately transparent but had a 
greater mismatch of refractive indices, exhibited much lower 
transmission ratios of 0,10 and 0,04- respectively, for 
fibres perpendicular and parallel to slit. This shows that 
a visual inspection for transparency is inadequate. A 
scientific method, such as the one described above, should 
be used for precise and quantitative measurements of 
transparency of materials* 

3, 4-, 2 Resolution of Lines Through Model Materials 

It was pointed out earlier that the sharp boun- 
daries of an object may appear smeared when viewed through a 
composite model whose constituents do not have well matched 
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refractive indices. This effect has been examined by viewing 
closely spaced straight and parallel lines through the models. 
Two sets of equispaced lines were drawn inside a circle with 
lines m one set being perpendicular to -chose m other set 
(Figure 3.5 a). These lines were viewed through the circular 
discs of unidirectional and bidirectional model materials 
(Figure 3*5 b and c). Direction of fibres has been indica- 
ted by arrows. It is clear from the figures that the lines 
are easily distinguishable from one another m either case. 
Another unidirectional material with a small mismatch of 
refractive indices shows a different behaviour (Figure 3*5 d) . 
The lines which are perpendicular to the fibres still do not 
appear smeared while the lines parallel to fibres are affec- 
ted to an extent that they can not be distinguished from one 
another. 

For a quantitative measurement of the above effect, 
the parallel lines, with distance between two consecutive 
lines changing according to a logarithmic scale, \/ere 
viewed through unidirectional models. The spacing between 
the lines varied from 0.1 mm (pitch 0,2 mm) to 1,0 mm 
(pitch 2,0 mm). The behaviour of model material has been 
compared with the other inferior material. For the supe- 
rior model material the lines could be clearly distinguished 
from each other by nacked eye even for smallest spacing. 
Moreover, the observation is independent of orientation of 
fibres with respect to lines. A slightly inferior material 
exhibited a limit of resolution of lines which changed with 



Figuri 3.5 View of parallel lines through (a) no;model,(b) 
unidirectional model material, (c) bidirectional 
model material and (d) a material with small 
mismatch of refractive indices. Arrows indicate 
the direction of fibres. 
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the angle between the fibres and the lines. The variation 
m the limit of resolution of linos has been shown as a 
function of fibre orientation in Figure 3.6. The highest 
resolution of about 4 lines per mm (spacing 0,125 mm) is 
obtained when fibres are perpendicular to the lines. The 
limit of resolution markedly leducos to only one line per mm 
(spacing 0.5 mm) for fibres parallel to lines. The diffe- 
rence in resolution of lines through the two materials is 
due to optical irihomogeneity of the second material winch 
arises due to the difference m refractive indices of its 
constituent materials. 

3.4,3 Image Forming Characteristics 

Optical inhomogeneity m the path of light causes 
its deflection and consequently, a directed beam of light 
gets deviated. As discussed earlier, optical homogeneity 
of the model material, which is directly related to the 
accuracy of matching of refractive indices of its consti- 
tuents, can be examined by studying the image forming 
characteristics of the material. For this purpose the 
optical bench, shown schematically m Figure 3*2 and dis- 
cussed earlier, has been used. With this arrangement of 
optical elements a sharp image of the slit is formed on 
the screen (Figure 3*3 4). The influence of placing the 
model between the second lens and screen is examined heie. 
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Photographs m Figure 3.7 show images of the light souice 
with models placed m position indicated in Figure 3*2. 

The images m Figure 3.7 a, b and c correspond to the uni- 
directional model materiel vitn fibres oriented at 0, k-5 and 
90 degree to the slit. The image m Figure 3.7 d corresponds 
to the bidirectional model material. It can be clearly seen 
that there is a negligible distortion of images when the 
superior model materials are placed in the path of light. 

The main image can be clearly distinguished from the diffused 
light which is observed in the direction perpendicular to 
fibres. This shows that the refractive indices of the glass 
fibres and the polyester matrix have been matched quite ade- 
quately and the materials may be regarded as optica] ly homo- 
geneous . 

The effect of mismatch of refractive indices can 
be seen m Figure 3*8. Photographs m Figure 3.8 a and b 
correspond to a glass-polyester composite with co/nposition 
of polyester only slightly different from that of superior 
model material. Images m Figure 3.8 c and d correspond 
to a glass-epoxy composite which has still larger mismatch 
of the refractive indices. It rnay be pointed out that on 
visual inspection all three materials look adequately 
transparent and may be thought of as suitable model mate- 
rials for photoelastic applications. However, capability 
of retaining a directed light beam m its original form 
makes a material superior than the others. 
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and (c) 90® and (d) image through bidirectional 
model material. 
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Brie effect of mismatch of refractive indices 
becomes more acute when the incidence of light on the mate- 
rial surface occurs at an oblique angle* It has been 
observed that the image of the straight slit becomes curved 
by oblique incidence on models with mismatched refractive 
indices of the constituents. However, no distortion of the 
image observed by oblique incidence on the superior model 
materials . 

3.4.4 Photoelastic Effects 

The model materials nave been examined for photo- 
elastic effect by placing their circular disc models m a 
circular polanscope employing a diffused sodium light 
source. The discs were loaded m diametral compression 
with fibres oriented at 45° to tile loading axis. The 
resulting iso chromatic fringe patterns for unidirectional 
and bidirectional models have been shown m Figure 3*9 a 
and b respectively. The isochromatic fringe patterns for 
another glass-polyester composite model with slight mismatch 
of refractive indices and a glass-epoxy composite model 
with larger mismatch have been shown m Figure 3*9 c and d 
respectively. It can be clearly seen that the fringes for 
the superior model materials (Figure 3.9 a and b) are much 
better defined than for other materials (Figure 3*9 c and d) . 
Some nonuniformity of fringes parallel to the fibres occurs 
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model (b) bidirectional model (c) glass- polyester 
model with small mismatch and (d) glass-epoxy 
model with large mismatch of refractive indices. 
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due to nonuniform distribution of fibres so that photoelastic 
sensitivity of the material changes from point to point. 

From the fringe pattern shown m Figure 3*9 c, it may appear 
than the inferior glass-polyester composite is also an ade- 
quate photoelastic material but due to undesirable charac- 
teristics (Figure 3*5 d, 3*8 a and b) as o result of larger 
mismatch of refractive indices, this material may not be 
suitable for more sophisticated analyses, 

3.5 MECHANICAL AND PHOTO ELASTIC CHARACTERISATION 

Mechanical and photoelastic properties of bire- 
fnngent materials play an important role in their selection 
as model materials. The mechanical properties of orthotropic 
materials, which are relevant to photoelastic analyses, are 
the elastic moduli along the two principal material axes, 
shear modulus and the major Poisson's ratio. These proper- 
ties were evaluated m uniaxial tension tests on specimens 
with different fibre orientations. The two elastic moduli 
and Poisson's ratio are directly determined from tests on 
specimens with fibres oriented at 0 and 90 degree to the 
load axis. The shear modulus has been indirectly deter- 
mined by measuring off axis elastic modulus on a 45 degree 
specimen and using transformation equations given m 
Appendix A. All the tests were performed on an Instron 
machine where the load was measured by a load cell and 
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strains were measured by mounting electrical resistance 
strain gauges on the specimens* The mechanical properties 
of the unidirectional and bidirectional materials are given 
m Table 3*1* It may be observed that the ratio of the two 
elastic moduli is 2.7 for unidirectional material which indica- 
tes that the material possesses adequate mechanical aniso- 
tropy. Anisotropy of the material may be changed by chan- 
ging volume fraction of fibres. 

The property that determines the photoelastic sensi- 
tivity of a model material is the stress -f ringe value or 
equivalently, the strain-fringe value. For an isotropic 
material stress-fringe value f and strain fringe value 
f £ are independent of direction and are defined as 


a <5 



where 

0 p - = principal stress difference 

= principal strain difference 
N = relative retardation or fringe order 

For orthotropic materials, the f a an d f e change with 
the direction. The theory of photoelasticity for ortho- 
tropic material will be discussed m the next two chapters 
but it may be mentioned here that photoelastic characteri- 
zation of an orfchotropic material requires three independent 


photoelastic eonstants that is either the stress-fringe 
values or the strain -fringe values defined as follows 
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where subscripts L and T refer to the longitudinal and 
transverse directions m the manorial. The stress-fringe 
values, % L and f a rnj are determined by applying uniaxial 
stresses in the L and T directions respec cively. The strain 
fringe values f e L and f £ ^ mav be determined by applying 
uniaxial strains m the L and T directions respectively. 

Hie and fg^rp are determined by applying the shear 

stress and shear strain respectively m the L and T direc- 
tions. In the present studies the f 0 j and f 01 p were deter- 
mined directly using tension specimens with fibres oriented 
at 0 and 90 degrees to the load axis. The strain-fringe 
values f £ y and f £ p are difficult to determine directly due 
to the experimental difficulties m applying the uniaxial 
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strains. Therefore, strain-fringe values f ^ and f ^ 

and the stress fringe value were evaluated indirectly 

using theory of photo-orthotropic elasticity to be discussed 
m the following chapters. The values of the six photo- 
elastic constants for the model materials axe given m Table 
3.2. These values are quite comparable with the photoelastic 
constants of other orthotropic materials reported m 
literature jl9, 2Q, 25 
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Table 3*1 Elastic Properties of Model Materials 



i , 

| e t 

! Q ! 

1 /LT { 

1 6 , 2 • 

’ 10 U/m 1 

V LT 

Material 

> 6 ? 
110° rl/rri 

r f, o 

! 10 N/m 
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Bidirectional 
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7734.5 
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Table 3.2 Photoelastic Properties of Model Materials 
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Bidirectional 

90.5 

80.5 

35.8 

13.634 

12.624 

10.137 


* Wavelength of light, X = 5893 A 



CHAPTER 4 


EXACT THEORIES OF PHOTO-ORTHCTROPIC-SLASTICITY 

4.1 INTRODUCTION 

The second requirement for the application of 
transmission photoelasticity to the stress analysis of 
composite materials is to develop an adequate theory 
through which photoelastic observations can be interpreted 
and translated into desired stresses or strains. Micro- 
scopically, glass fibre reinforced composites are hetero- 
geneous with each constituent being homogeneous and 
isotropic. On a macroscopic scale, the fibrous composites 
are often assumed to be homogeneous but anisotropic. There- 
fore, the photoelastic theories for composite materials 
have been developed by considering their behaviour either 
on microscopic or macroscopic scale. In the former case, 
the photoelastic response of each constituent of composite 
material is taken into consideration to predict the overall 
response of the composite. In the latter case, the behaviour 
n£ oxjmpo-site material is assumed to be analogous to aniso- 
tropic crystals and therefore, their photoelastic response 
is predicted using photoelastic theories of crystals. In 
this chapter some of the existing theories of photo - 
orthotropic-elasticity will be reviewed first. Then, an 
exact strain-optic law will be derived with a view that it 
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may be later modified to remove limitations of the existing 
theories of photo-orthotropic-elasticity. 

4.2 ASSESSMENT OF EXISTING THEORIES OF 
PHOTO-ORTHOTROPIC-ELASTICITY 


Pih and Knight [17] were probably the first to 
propose a stress-optic law for orthotropic materials. Their 
approach was based on the distribution of stresses between 
the constituents of a composite. However, their theoretical 
predictions and experimental results did not show good 
agreement. Therefore, their approach was not pursued further 
by other investigators. Sampson 1 1 9 3 developed the concept 
of a Mohr circle of birefringence by simple analogy. He 
postulated that the orthotropic materials can be characterised 
with three principal stress-fringe values f 0L , f 0 m} f 0LT j 
and that the stress system and produces three 

components of birefringence as follows: 
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where the subscripts L and T refer to the principal material 
axes (for example, the directions parallel and perpendicular 
to fibres in a unidirectional composite). The resultant 
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birefringence, then, is: 
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The principal stress-fringe values f CT -, f arp and f a ^ can 


be experimentally determined through three independent tests 
to be discussed later in this chapter. 

An off-axis stress-fringe value, f a0 , under -uni- 
axial state of stress, nay be defined as s 
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where is the uniaxiam stress applied at an angle, 9, to 
the longitudinal direction. Using the laws of stress trans- 
formation, f _ may be obtained in terms of f T , f m and 
’ o y o ij o i 

ii as * 

2 

o f _ _ p 2 f_ - p 

f aQ = f aL C ( cos 9 - ~ sin 9 ) + — ^ sin 2 9] 

° T f aLT 

C^+»^+ ) 

The optical isoclinic parameter, 0 (that is, angle that an 
isoclinic fringe makes with the fibre direction) can be 
computed as a function of principal stress direction, 9, for 
uniaxial stress % 

f 0T sin 2 9 

tan 2 0 = j 2 Of. 5) 

f crLT ^ cos e “ f al/ f aT • sin ©) 

Sampson [19 1 supported his theory by limited experimental 

results which were obtained by testing unidirectional 
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composites under uniaxial state of stress. 

Dally and Prabhakaran [20 ] , like Sampson [19 ], 
assumed that the photoelastic response of orthotropic mate- 
rials can be predicted through chree independent stress- 
fringe values. They derived a stress-optic law for uniaxial 
state of stress. Conceptually, their approach is b r sed on 
the finding out the stresses, and hence the contributing 
birefringence in each constituent. The net birefringence is 
an algebraic sum of the contributions bjr each constituent. 

The accuracy of results is greatly influenced by the accuracy 
with which the stresses in constituents can be evaluated. 

This becomes a limitation of the approach particularly because 
it is very difficult to predict transverse stresses accurately. 
Moreover, they did not discuss the interpretation of isocli- 
nics in orthotropic birefringent materials. Experimental 
and theoretical results of Dally and Prabhakaran closely 
agreed with the stress-optic law formulated by Sampson. 

Bert [22] modelled fibrous composites by orthor- 
hombic crystalline system and applied the theory of photo- 
elasticity for anisotropic crystals which was initially 
propounded by Pockels [231 and later modified by Bhagavantam 
[2?+]. According to this theory the symmetric second order 
birefringence tensor is related to stress or strain tensor 
through a fourth order photoelastic-property tensor. Using 
this theory Bert derived a stress-optic law for orthotropic 
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materials subjected to plane stress and sieved that the 
concept cf Mohr circle of birefringence, as proposed intui- 
tively by Sampson [19 3, is a direct result of tenscrial 
nature of birefringence. The interpretation that the iso- 
clinic angle represents principal birefringence direction, 
is also consistent. More recently Knight and Pint 2 >3 using 
Bhagavantam-Pockels theory derived a general stress-optic 
law and then simplified it to a two dimensional case. Their 
formulation is general and accounts also for initial bire- 
fringence which is usually present in two-phase materials. 
However, their simplified equations are of the same form as 
those of Sampscnf 19 3 and Bert [223. 

It is clear from the above discussion that most 
of the existing theories of photo-ortho tropic-elasticity lead 
to stress-optic law of the form proposed by Sampson. Some 
experimental results reported in the literature indicate 
that the photoelastic response of orthotropic materials can 
be accurately predicted through three independent photo- 
elastic constants. It was, however, decided to assess this 
form of stress-optic law through the experiments conducted 
on superior orthotropic model materials whose development 
has been discussed in the previous chapter. 

Photoelastic response of model materials was stu- 
died on tensile specimens of unidirectionally and bidirec- 
tionally reinforced materials. The fibre orientations were 


riept at, 0, 1 5 , 30, 45, 60, 75 £.na 90 degree to fee longitudi- 
nal eag' of the specimens, "lie steelmens were gradually 
1 Octde a in uniaxial tension anc. ocserved in a transmission 
poiariscopu for isoclinic and isochrcmatic measurements. 

The resulting isoclinic parameter and i sc chroma tic fringe 
order was noted down ..or ea. i inc-'ement cf load. The stress- 
fringe value of a specimen was obtained by measuring slope 
of the plot of the fringe order (IT) as a function of applied 
stress. A typical fringe order vs stress curve is shown in 
Figure 4.1 for the specimen with fibre orientation cf 30 
degree. The value of f a ^ T was determined using an indirect 
procedure in which an off-axis stress-fringe value is 
measured and the value of f is calculated from So. (4.4). 

alii 

In general any off axis specimen can be chosen for the pur- 
pose but the effect of errors, which might have been made 
in the measurements of f gL and f can be minimised, if not 
completely eliminated, through proper selection of off-axis 
specimens. It can be seen from the Eq. (4.4) that f aLrp will 
be related to only f c g provided the term in parenthesis is 
zero. When this is achieved f 0 - ^ becomes s 


f aLT 



sin 2 9 q 


wher tan^ 9 = —14 (4,6) 

° 1 a L 

Angles 0 O for the unidirectional and bidirectional mate- 
rials were found to be 33 and 43 degrees respectively. 




in the eresent case f T -t was, therefore, calculated from 

aLi 7 7 

a 30 degree specimen for unidirectional material and from 
45 degree specimen for bidirectional material. 

The stress-fringe value has teen plotted as a 
function of fibre orientation, 9, in Figure 4.2 for both 
model materials. The difference observed in the values of 
f a _ and f aT for bidirectional material is expected because 
the amount of reinforcement in the two directions is not 
equal. Also shown in Figure 4.2 are the theoretical pre- 
dictions based upon Eq. (4.4). It is very clear that the 
experimental results agree very well with the theoretical 
predictions. Isoclinic angle, 0^, as measured from the 
direction of load, has beer, plotted as a function of fibre 
orientation, 9, in Figures 4.3 and 4.4 for unidirectional 
and bidirectional model materials respectively. Isoclinic 
angles predicted by Eq. (4.5) have also been plotted. It 
may be mentioned that Eq. (4.5) gives angle from the fibre 
direction where as the angles plotted in Figures 4.3 and 4.4 
have the load direction as the reference. Also shown in 
the figure are the directions of principal stresses and 
principal strains. It is observed that the experimentally 
measured values of isoclinic angles are in good agreement 
with the theoretical predictions. Further, the isoclinics 
do not coincide with either the directions of principal 
stresses or strains which means that the principal 




Fiber-orientation, ^(degree) 

Figure 4.3 Isoclinic angle as a function of the fib 

for uniaxial loading of tensile specimen 
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"birefringence directions in orthotropic material are, in 
general, not the same as the principal stress or strain 
directions. However, the isoclinics are closer to the 
directions of principal strains than the principal stresses. 
Also, the directions of principal strains show the similar 
variation as the isoclinics for both the materials. It may- 
be mentioned that the initial birefringence in orthotropic 
materials will influence the isoclinic angle. In the present 
case, however, the initial birefringence was very small In 
magnitude. And its influence was further minimised by 
measuring the isoclinic angle at higher stresses when the 
isoclinic angle becomes practically constant. 

It is now quite clear from the above discussion 
that photoelastic response of orthotropic materials can be 
accurately predicted through three independent photoelastic 
constants. It is, however, not possible to evaluate, 
through the existing stress-optic laws, the individual 
stress or strain components or the difference between their 
principal values from the birefringence measured in a polari- 
scope. Similarly it is possible to predict the optical 
isoclinic angle for a known state of stress but it is not 
possible to predict the direction of principal stress or 
strain from the measured optical isoclinic angle. Thus 
the existing theories of photo-orthotropic-elasticity are 
not capable of translating photoelastic measurements into 
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desired stresses or strains. This has been the aim of present 
investigations to develop a theory of photo-ortho tropic- 
elasticity which will overcome the above limitations in the 
applications of transmission photoelasticity to the problems 
of stress analysis of ortho tropic materials. 

Some approximate laws have been proposed by inves- 
tigators to overcome the limitations of theories of photo- 
orthotropic-elasticity discussed above. Pipes and Rose t 26] 
proposed the use of strain-optic law in place of stress- 
optic law. They assumed that isotropic strain-optic law 
(with a single strain-fringe value) is applicable for low 
fibre-volume-fraction composites where the matrix dominates 
the optical response of the material. They reasonably 
verified this assumption for composite with low fibre-volume- 
fraction. However, usefulness of their work is limited 
because the real composites often have high fibre volume 
content and because they did not indicate the magnitude 
of error that may arise due to the simplifying assumptions. 
Their work did show a new direction in which efforts could 
be made to make transmission photoelasticity more useful. 
Prabhakaran t 27 1 proposed a strain-optic law with three 
principal strain-fringe values and hence it has the same 
drawbacks as the stress-optic laws discussed earlier. In 
the same paper Prabhakaran has proposed a simplified strain- 
optic law with a single strain-fringe value which is taken 
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to be the arithmatie average of the three principal strain- 
fringe values. The approaches of Pipes and Rose and of 
Prabhakaran are really some kind of ad hoc approaches in 
which the strain-fringe values are calculated from the 
stress-fringe values through the composite stress-strain 
relations. The derivations have neither mathematical rigour 
nor a fundamental basis of the material behaviour. Their 
formulations failed to provide a method of taking the errors 
into account due to approximations involved. 

In the present investigations it was decided to 
develop an approximate strain-optic law. The motivation 
behind picking jap the strain-optic law and not the stress- 
optic law for approximations \^as that in orthotropic mate- 
rials, the optical isoclinic angles are closer to the 
directions of principal strains than the principal stresses. 
Secondly there are some indications that the ortho tropic 
materials exhibit lower degree of orthotropy with respect 
to strain-optic behaviour than that with respect to stress- 
optic behaviour. The approximate strain-optic law was 
derived from an exact strain- optic law which in turn was 
derived from the Bhagavantam - Pockels theory of photo- 
elasticity of crystals. The approximations have been made 
in a manner which minimises the error. Moreover, a method 
has been developed to estimate the error with a reasonable 
degree of accuracy. 



83 


If. 3 EXACT STRAIN-OPTIC LAW 


It is generally agreed that macroscopic behaviour 
of fibrous composites can be accurately predicted using 
theory of anisotropic elasticity. It, therefore, seems 
appropriate to assume' that the photoelastic behaviour of 
composites is analogous to that of anisotropic crystals. 

The theory of photoelasticity of crystals was first proposed 
by Pockets 1231. Later on Bhagavantam E 4l ] corrected cer- 
tain errors in it. In the established theory of crystalline 
photoelasticity [24] , birefringence is shown to be a 
symmetric second order tensor like stress or strain. The 
birefringence tensor N. . is related to the strain tensor 

% d 

( e^) through the following relation? 


N i3 ~ P ijkl E kl 


(i, j = 1, 2, 3) 


(4 .7) 


where is a fourth order strain-optic property 

tensor like the elastic coefficient tensor. However, 
unlike the elastic coefficient tensor, the strain-optic 
property tensor is generally unsymmetric. The tensor, 
p ijkl’ in § eneral , consists of 36 independent strain- • 
optic coefficients for a three dimensional anisotropic 
birefringent material. 

For simplicity Eq. (4.7) may be written in the 
contracted notations as: 


N„ 


P e 
qr r 


(q, r = 1, 2, ... 6) (4.8) 
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Behaviour of an ortho tropic material, such as 
fibrous composite, can bo considered analogous to that of 
an orthorhombic crystal having three principal material 
symmetry axes. Bor orthorhombic system the strain-optic 
relation. So. (4.8) takes the following form when written 
in the matrix notation s 


(4.9) 


It may be noted that there are 12 independent strain-optic 
coefficients for an orthotropic material), compared to only 
9 elastic constants. 

Using the same reasoning as used for the analysis 
of mechanical behaviour of fibrous composites l 35^, it can 
be assumed that the unidirectional material, on a macrosco- 
pic scale, is transversely isotropic, with the plane of 
isotropy being noimal to the fibres. If the fibre direc- 
tion is taken as axis 1, Eo. (4.9) becomes: 

’’Here onwards e j +5 and. e g will be denoted in terms 
of shear strains as 1 23’ Y 31 andY -i2 * 




8 5 


r* 

— 





~ I i 


u i 


P 11 

P 12 

P 12 

0 

0 

0 . 


e i 

H 2 


P 21 

P 22 

P 23 

0 

0 

0 


E 2 

H 3 


P 21 

P 23 

P 22 

0 

0 

0 


e 3 

\ 

M 

0 

0 

0 

P 22“ P 23 

“2 

• 0 

0 

1 


Y 23 

%: 
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0 

0 

0 

P 66 

o : 


Y 13 

h 6 

1 

0 

0 

0 

0 

0 

p 66 


Y 12 








1 

~ — 


(4.10) 


A large number of real -life problems are idealised 
to be plane stress problems. The strain-optic relation for 
plane-stress becomes s 


N 1 

■ 

P 11 

P 1 2 

P 1 2 

0 


£ i 

k 2 


P 21 

p 22 

p 23 

0 


e 

2 

H 3 

i 


P 21 

P 23 

P 22 

0 


e 3 

N 6 

j 

0 

0 

0 

p 66 


j r 12 

^ 


— 



— 



It may be noted that, for a plane stress case, the out-of- 
plane strain, e is, in general, not zero and will affect 
the in-plane birefringence. Hie ref ore, the terms in the 
third row and third column of the matrix in Eq. (4,10) 
can not be dropped. However, for a plane stress case of 
transversely isotropic composite, the strain, normal 

to the plane of composite lamina can be written in terms 
of in-plane strains, and, as follows; 
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C-v 12 (1 + v 23 ) £l 4 - ( v 23 + v 12 v 21 )c 2 ] 


12 21 


(4-.12) 


or 


e 3 “ 


K £ - K„ e„ 

1 1 2 2 


(4.13) 


where 


K i « 


K 2 = 


V 12 (1 + V 2 3 } 
1 - v 12 v 21 

V , V V 
23 12 21 
1 - v 12 v 21 


(4.14) 


The birefringence components in the plane 1-2 can be 
written in terms of strain components from Eq. (4.11) as 
follows i 


N„ - N, 


(P-J-^ — PoJ ^ ~ (Poo ~ PlO^O + (P-,0 “ PoQ^ 


21 


22 


12 ' 2 


N 6 = P 66 Y 12 


12 ^ 23 ' 

(4.15) 


Now, by substituting la. (4.13) into Eqs. (4.15), the in- 
plane birefringence components can be written in terms of 
the in-plane strains only as: 


N 1 - N 2 = b 1 e 1 - b 2 e 2 


% = h Y 12 


where 


_ (P 11 - P 21 ) - (P 12 - P 23 ) 
= (P 22 - P^ 0 ) + K 0 (P 19 - P 0 ^) 


12 


12 


23' 


66 


(4.16) 


(4.17) 
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The quantities observed in a polariscope are the resultant 
birefringence and the direction of principal birefringence. 
The resultant birefringence is the difference between the 
principal components of birefringence. The resultant bire- 
fringence represents relative retardation and is measured 
through the isochromatics which are lines of constant 
relative retardation. The direction of principal birefrin- 
gence is represented through the angle it makes with a 
fixed reference axis. The angle is called isoclinic angle 
and the lines of constant isoclinic angles are called iso- 
clinics. 

Birefringence components in an arbitrary direction 
can be obtained from known magnitude and direction of 
principal birefringence components through the second-order 
tensor transformations. It can be easily shown that the 
brief ringence components given in Eq. (4.16) can be written 
in terms of resultant birefringence, N (the difference of 
principal birefringence components N , and N^,), an( 3- the 
isoclinic angle, 0, measured from the fibre direction, as 
follows ; 

N 1 - N 2 = N cos 2 0 = b 1 c 1 - b 2 c 2 

2 N 6 = I sin 2 0 = 2 b 3 Y 12 (4.18) 

The desired strain-optic law can now be written 
by separating the resultant birefringence and isoclirdc 
angle : 
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where, f^.j, is known as the transverse normal strain- 
fringe value* 

(iii) b^ can he determined when only pure shear strain, 

Y lt , is acting on the material. Under these 
conditions : 



_J _ N 

2 f £LT Y LT 


(4-. 23) 


With the definitions of strain-fringe values given 
by Eqs. (4-. 21), (4-. 22) and (4-, 23), the strain-optic law 
can be written in a more familiar form as : 


N 


l ( 




Cm 2 Y -r rp 2 1/2 

r 1 ) + ( ) ] 

l£r l £ LT 


C4-.24-) 


and the isoclinic parameter, 0, becomes: 


0 =4 tan" 1 


f eLT 


(4-. 25) 


•EL 


'£T 


The above form of strain-optic law is convenient when 
strains acting along material symmetry axes are under con 
siderations. In many stress analysis problems, it may be 
necessary to consider the principal strains and their 
orientations. If the principal strains are denoted by 
and where the direction of, makes an angle, ©, 
with the fibre direction, the strains along L and T 
directions can be obtained by transformation as : 
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2 . 2 

e r = e cos 0 + e sin 0 

L p q 

. 2 2 
e T = e p sin 9 + Eq cos 0 

Y lt = ( e - £q) sin 20 


Of. 26) 


Eqs. (4-. 26) and (4.24) can be combined to obtain a general 
strain-optic law for biaxial strain- field as s 


N 


e r>, p s ct p p e q p 

•p — { t (cos^ © + ~ sin^ 0) - ( sin^ 0+7— cos 0) 

EL n D 


SL 


2 e f 2 1/2 

J + [ (1 - 7^) b ' £L sin 2 0] > 


%' f £LT 


(4.27) 


The expression for isoclinic parameter, 0, becomes; 


(1 - 


0 = ^ tan”" 5 


cos 2 0 


_EL 


sin 2 0 
- 

eLT 


. 2 e_ 
sm 0 4 , 

~T~ + e { 


sin 2 0 


2 

cos 0 


f f r e v f f 

£L eT p £L eT 


) 


(4.28) 

It may be pointed out that the strain-optic law proposed 
intuitively by Prabhakaran t2 7] is in agreement with the 
exact strain-optic law given by Eqs. (4.24) or (4.27)* 
However, it was not clear in Prabhakaran 1 s approach if and 
how the effect of out-of-plane normal strain, e ^, is 
taken into account. The exact strain-optic law is also 
of the same form as the stress-optic laws formulated by 
Sampson C19], Bert [22] and Knight and Pih [2 51 all of 
which also involve three stress-fringe values. 
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It will be appropriate to simplify the strain- 
optic law r for the following two special cases of practical 
significance* 


(a) Strain-optic law for uniaxial stress 

For uniaxial state of stress the ratio of principal 
strains is given by the Poisson's ratio ass 



v 

pq 


(4.29) 


An off-axis strain fringe value, f £ Q, may be defined ass 


I ee “ N ~ 

The following expressions for, f 
by combining Eqs. (4.27), (4.28) 


S (1 + V Ba ) 

N 


(4.30) 


£ g, and 0, may be obtained 
, (4.29) and (4.30) s 


f <=e = C1 +V pq > W ' (cos ‘ e - V pq sln2 


. 2 
sin 

© - 

V 

cos 2 e) Jh 

] +[ (1 + v ) 



pq 


pq 




2 -1/2 


£L 

sin 

2 © 

i > 

(4.31) 

eLT 





-1 


(i 

+ V) 

1 e LT 

_© 


2 2 2 2 
cos 0 sin Q v / sin Q cos Q \ 

f eL f el pq f eL f eT 


(4.32) 
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(b) Strain-optic law for uniaxial strain 

In the case of uniaxial strain (i.e., e = 0), an 
off-axis strain-fringe value f| Q may be defined as; 


f £ ’g 


_p 

Li 


(4. 33) 


The following expressions for, f^ G , and 0, may be obtained 
by combining Eqs. (4.27) , (4,28) and (4.33) and substitu- 
ting e ^ = 0 ; 


ft 

e 9 


P = 


1 

2 


2 Q 

[ ( 40 1__£ . 
% L 


tan* 


cos 


sin 2 9 2 

f ' 

X £rp 

sin 2 Q 
7 ” 

eLT 

Q sin 2 0 


+ ( 


sin 2 0 
f eLT 

(4.34) 

(4.35) 


f eL f eT 

It may be pointed out that although the expression (Eq. 4.34) 


for off-axis strain-fringe value is similar to that of 
off-axis stress-fringe value (Eq. 4.4), the direct determi- 
nation of, j^ Q , is considerably more difficult than that of 

f In order to obtain a uniaxial strain-field the material 
a © 

has to be subjected to a biaxial stress-field. Moreover, 
the ratio of stresses in the two directions has to be 
adjusted precisely to compensate for the Poisson's effect. 
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4-. 3. 1 Experimental Determination of Principal 
Strain-Fringe Values 


Application of the strain-optic law derived in the 
previous section requires three principal strain-fringe 
values which have to be determined experimentally. As 
discussed earlier, their direct determination will require 
the application of biaxial stresses on the specimen. It 
was therefore decided to follow an indirect approach for 
their determination by applying uniaxial stress only. 
Following are the steps followed in the determination of 
the strain-fringe values. 


(i) A uniaxial stress, along fibre direction produces 
the principal strains along directions (1 - 2) as 
follows ; 


£ 1 = £ L 


e _ «_ v e 

2 ~ LT L 


( 4 . 36 ) 


The isochroma tic and isoclinic response can be obtained 
from Eqs. (4-. 19) and (4.20 ) j 


N = (b, + 4 2 V lt )S l 

0 = 0 


( 4 . 37 ) 


The strain can be calculated from the known longitudi- 
nal stress and elastic properties or may be directly 
measured by fixing electrical resistance strain 
gauges on the specimen. A plot of isochromatie fringe 
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order as a function of longitudinal strain is shown 
in Figure 4. 5. The straight line thus obtained has 
a slope, S.j, given by; 

s i = h + '°2 \i c^.38) 


(ii) A uniazial stress, a T , in the transverse direction 
will produce the following principal strains: 


e 1 ~ " V TL e T 


e 


2 


e 


T 


(4-. 39) 


The isochromatic and isoclinic response can be obtained 
from Eqs. (4.19) and (4.20); 

N = (b 1 V TL + b 2 ) e T 

0 = tt/2 (4 .40) 


A plot of isochromatic fringe order as a function of 
is shown in Figure 4-. 6 . The straight line thus 
obtained has a slope, S 2 , given by; 

S 2 = b 1 V TL + h 2 (4-. 4-1) 


The Eqs. (4-. 38 ) and (4-. 4-1) can be solved for b^ and b 2 as: 


5 1 - S 2 \T 
1 " V LT V TL 

5 2 - S 1 \L 

1 - v v 

1 LT TL 


(4.4-2) 
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The principal strain-fringe values are given in terms 


of b-j and b£ 
the, f £ ^ and 

by Eqa. (4,21) and (4.22). Therefore, 
f e j become : 

-p _ 

1 „ V V 

LT TL 

sL “ 

■p _ 

S 1 " S 2 V LT 

(4.43) 

1 ~ V LT V TL 

f eT ~ 

S 2 “ S 1 V TL 


(iii) The direct determination of, b^) or the shear strain- 
fringe value, f e will require application of pure 
shear strain or equivalently pure shear stress. In 
an Indirect method if the off-axis strain-fringe 
value, f^ q, can be determined, the f eLT , can be 
calculated from Eq. C 4- . 3^ ) • However, due to the 
difficulties in applying a uniaxial strain, it was 
decided to calculate the shear strain-fringe value 
from the measured shear stress-fringe value. In 
view of definitions of f e ^ and f^^ and the relation 
between shear stress and strain, the relation between 
f eLT and f aLT can s - ll0wn to ^e : 


cfLT 


eLT 


2 G. 


LT 


( 4 . 44 ) 


The method of obtaining Has already been dis- 

cussed. The constant, b^j is the reciprocal of, 
f as given by Eq. (4.23). So that! 


b 


3 


2 g lt 
f cLT 


( 4 . 45 ) 
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The relations between the normal strain-fringe 
values -.nd the normal stres '-fringe values can also be 
obtained as % 


e L 


1 - v v 

LT TL 

e l 

f V T.T 
o T h 


a L 
1 - v 


LT V TL 


e T 


E, 




(4.46) 


TL 


cr T a L 

The strain-fringe values and stress-fringe values 
for the superior, unidirectional model material have been 
given in Table 4.1. Some important observations on them 
as well as on the photoelastic constants reported by other 
investigators will be made in the next section. 


4.4. SOME IMPORTANT OBSERVATIONS 


In view of the results and discussion presented in 
the previous section the following important observations 
can be made i 

(i) The variation in the strain-fringe values, f e q, for 
a uniaxial stress and, f f £ Q, for a uniaxial strain 
as shown in Figure 4.7 is less than + 10 percent of 
their respective average values. This variation is 
much smaller compared to that in the stress-fringe 
value f a q, (Figure 4.2), which is nearly 100 per 
cent for a particular fibre orientation. 



Table 4,1 Comparison of Photoelastic Properties and Anisotropic Ratios 
of Different Model Materials 
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(ii) The principal strain-fringe values and stress- 

fringe values and their ratios have teen given in 
Table 4.1 for the superior unidirectional model 
material as well as for some other materials repor- 
ted in the literature. The strain-fringe values 
for the materials of other investigators have been 
calculated from the reported stress-fringe values 
and elastic constants using Eqs. (4.44-) and (4.4-6). 
The strain-fringe values directly reported by the 
investigators may not be in agreement with the 
values calculated in this manner because of the 
difference in strain-optic laws. The ratios f el/ f eT 
and f eT /f are indicative of the anisotropy with 
respect to strain-optic behaviour of materials, 
whereas the ratios are indica- 

tive of material anisotropy with respect to stress- 
optic behaviour. It is observed that the ratios 

f £ l/f gT and ^eT^eLT are closer "1*° than the 
ratios f 0L /f aT and f^/f^j. This indicates that 
the degree of orthotropy of materials with respect 
to strain-optic properties is considerably smaller 
than that with respect to stress-optic properties. 

(iii) In view of results shown in Figures 4.3 and 4.4, 
it is observed that the directions of principal 
birefringence (directions of isoclinics) are much 
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closer to the directions of principal strains than the prin- 
cipal stresses. The maximum angle "between the principal 
birefringence direction and the principal strain direc- 
tion is less than 5 degree where as the maximum angle 
between the principal birefringence direction and princi- 
pal stress direction is more than 15 degree. Moreover, 
the variation in principal strain direction with fibre 
orientation is qualitatively also similar to that in the 
isoclinic angle. 

The above observations strongly suggest that an 
approximate strain-optic law can be developed through which 
photoelastic observations can be interpretted without 
loosing much accuracy. An approximate strain-optic law 
has been developed in the next chapter. 



CHAPTER 5 


AN APPROXIMATE THEORY OF PHOTO-ORTHOTROPIC-ELASTICITY 
5.1 INTRODUCTION 

Exact theories of photo-orthotropic-elasticity 
have been discussed in Chapter 4. It has been established 
that the photoelastic response of orthotropic composites 
can be accurately predicted using exact stress-optic or 
strain-optic laws. The exact laws require three photo- 
elastic constants for the material which are determined 
experimentally. It was, however, pointed out that the 
exact laws are inadequate to interpret the photoelastic 
observations and translate them into the desired stresses 
or strains. An approximate strain-optic law will be 
derived in this chapter. A method of applying a correc- 
tion to the approximate strains will also be developed 
which reduces the error significantly. The reasons for 
picking up the strain-optic law and not the stress-optic 
law for approximations are based on the observations 
made in Chapter 4. It has been observed that the bire- 
fringent composite materials exhibit a smaller degree of 
orthotropy with respect to strain-optic behaviour than 
with respect to stress-optic behaviour. It has also been 
observed that the directions of optical isoclinics are 
much closer to the principal strain directions than the 
principal stress directions. 
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Before going into the process of approximation, 
it will he appropriate, at this point, to examine the 
concept of anisotropy. Anisotropy literally means devia- 
tion from isotropic behaviour. Isotropy and anisotropy 
are not general properties of a material but they refer 
to the material behaviour with respect to a specific pro- 
perty only. Anisotropy of a material with respect to a 
particular property does not necessarily imply that the 
material will be anisotropic to the same degree with 
respect to any other property. It is in this context that 
the photoelastic composites are observed to be anisotropic 
to a smaller degree with respect to strain-optic behaviour 
than with respect to stress-optic behaviour. The degree 
of anisotropy with respect to strain-optic behaviour, as 
observed in Ghapter 4, is small enough to justify represen- 
tation of strain-optic properties by an isotropic property 
tensor. 

5.2 APPROXIMATE STRAIN-OPTIC LAW 

To device an approximate strain-optic law it is 
assumed that the fourth-order strain-optic property 
tensor is isotropic having only two independent strain- 
optic coefficients. These two constants are analogous 
to the two independent elastic constants for mechanically 
isotropic materials. The independent photoelastic cons- 
tants are calculated from the constants of exact strain- 
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optic property tensor such that the error between the 
exact and approximate property tensor is minimum. This 
has been achieved by following the approach of Fedorov [ 38 ] 
used in the simplification of the study of stress waves 
propagation in crystals. 

The approximate strain-optic property tensor, 
P i;jkl’ assumec ^ 4o "b e the following most general fourth- 
order isotropic tensor? 


*lm = 0 5 i 3 5 ki + a ( hi +s ii V (5>1) 

where ^ is the Kronecker delta and a and c are 
the independent strain-optic coefficients. The square 
error, E, between the approximate property tensor given 
by Eq. (5*1) and the exact one given by Eq. (4.7) is? 


E ~ ( P ijkl “ P i3kl^ (P ijkl ” P iokl^ 


(5.2) 


The constants a and c are obtained by imposing the 
following conditions to minimise the error E: 


3E 

9c 


= 0 


9E 

3a 


= 0 


(5.3) 


Substitution of Eqs, (5.1) and (5.2) into Eq. (5-3) and 
solving for a and c yields ; 
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15 


(2 P - - 


H3 3 


a “ 30 (3 P iji3 


P i3'i3^ 


P ii3 3 ^ 


(5. 1 *) 


The coefficients P^^ in Eq. ( 5- 1 +- ) can he written in the 
following contracted notations consistent with those used 
in the Chapter for transversely isotropic composite as: 


c 

a 


T5 ( P 11 + p 22 + ^ P 12 + ^ P 21 + 5 p 23 ~ h P 66 5 


i r 


30 


2 (P^ - P 21 ) + 2 (P 22 - P 12 ) + 12 P 66 


+ ? ( p 22 * F 23 ^ ] 


(5.5) 


Approximate strain-optic relation for a plane stress case 
can be written in terms of a and c as follows : 


N 1 

L 

2a+c c c 0 

1 

e 

1 

U 2 


c 2a+c c 0 


e 

2 

k 3 


c c 2a+c 0 

! 

e 3 

HS 


0 0 0 a 

I 

Jl2__ 


(5.6) 


The in-plane birefringence components can now be obtained 
from Eq. (5.6) and written, as follows, in terms of the 
resultant birefringence and birefringence direction or 
isoclinic angle, 0, as was done in previous chapter: 



■l 

- H 2 

= N 

cos 2 0 

= 2 a 

( e. 


2 Kg 

II 

sin 2 0 

= 2 a 

Y 12 
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(5.7) 


The resultant birefringence and isoclinic angle can be 
separated as : 

2 2 1/2 

N = 2 a [ ( - ^ 2 ) + ( T^2^ 3 ( 5 * ^ ) 

Y 

tan 2 0 = — 7— (5.9) 

1 " 2 

It may be recognised that the square root term in the 
Eq. (5-8) represents the difference in the in-plane princi- 
pal strains. Representing the principal strain difference 
by ( £ p - e ^),the approximate strain-optic law may be 
written as : 


e - e = N f _ (5.10) 

p q e 

where f e is the strain-fringe value of the approximating 
body and is equal to l/(2a). The Eg. (5-9) implies that the 
directions of principal birefringence coincide with the 
directions of principal strains. In other words the iso- 
clinic angle is the same as the angle between principal 
strain and the reference axis. This is very significant 
observation and, along with the approximate strain-optic 
law (Eq. 5.10), this will significantly simplify the 
problem of obtaining individual strain components. 
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Hie approximate strain-optic law, given by 

Eq. (5»10)j contains only one photoelastic constant. Its 

value should be calculated from Bqs. (5-5) so that the 

error between approximate and exact theory is minimum. 

The Eqs. (5»5) gives the value 'a* (or equivalently 2 ~f~) 

in terms of six fundamental strain-optic constants namely 

P 11? P 12 , P 21 , P 2 2 ? P 23 and p 66* 13:10 strain-fringe value, 

f , can be calculated if the six constants are known or 
e 

if the values of their four groups as indicated in Eqs. 

(5.5) are known. Ey regrouping constants on the right 
hand side of Eqs. (4.17), it can be easily shown that the 
four groups of constants namely (P^ - P 2 -|)j ( p p2 ” p ']2' 3? 

P ^ and (P 22 - P 23^ are common to Eqs. (4.17) and Eqs. (5*5) • 
Thus Eqs. (4.17) will be quite helpful in calculating the 
values of these four groups. However, Eqs. (4.17) is set 
of only three Independent equations and can be used to 
calculate only three groups. Therefore, an additional equa- 
tion is needed to evaluate one of the groups or a combina- 
tion of them. To this end attention has been focussed on 
the group (P 22 “ ^23^* 

It follows from Eqs. (4.11) that the term 
(P 22 “ P 23> can be evaluated through an experiment in which 
the model is a thin 2-3 plane so that the birefringence 
(N 2 - N^) is measured and related to strains in the 2-3 
plane. However, a number of practical difficulties are 
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likely to be encountered, in carrying out such an experiment 
due to inhomogeneity of the material (glass fibres being 
discrete units) . It was, therefore, decided to evaluate 
the tern (P 22 - P^) through theoretical calculations. 

It is realised that the 2 - 3 plane behaves as an isotropic 
plane to in-plane loads so that a theoiy applicable to 
microscopically heterogeneous but macroscopically isotropic 
material can be applied to this case also. Richard and 
Young [39 1 using a self consistent approach, developed 
a theory for and demonstrated its applicability to an 
isotropic composite namely glass beeds reinforced polyes- 
ter resin with varying glass content. Therefore, it seems 
reasonable to use this theoiy for calculating (P 22 ~ P 2 ^). 
Following the theory, it can be shown (refer to Appendix C) 
that! 



E 

1 ( y 

1 + v 23 e 



+ x m 



am 


) 

(5.11) 


where V is glass content or volume fraction; f _ and 
g J a g 

^ m are the stress-fringe values for glass and matrix 
material respectively; X g and X m are respectively the 
ratios of principal stress difference in glass and matrix 
to that in the composite. These ratios X g and X^j can 
be evaluated in terms of elastic constants by following 
the approach of Richard and Young. The elastic constants 
of the isotropic composite can be determined experimentally 
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or estimated through self consistent approach adopted by 

Richard [ 4-0 ] . Expressions for Xg and and other 

details concerning derivation of Eq. (5.11) have been given 

in Appendix C. Now, that the group (P 22 " P 23^ is im ° wn 

the other groups required to calculate strain-fringe value, 

f , can be calculated using Eq. (4-. 17). In this manner 
£ 

the strain fringe value was found to be % 

-6 


9.90x10 m/fr 


( 5 . 12 ) 


5.3 ASSESSMENT OF STRAIN-OPTIC LAWS 

Photoelastic response of a material to any state 
of strain can now be predicted using exact or approximate 
strain-optic law. Predictions have been made for the 
superior unidirectional model material developed for photo- 
elastic applications and whose photoelastic properties are 
given in Table 4-.1. Wherever possible, calculations have 
been made for other materials reported in the literature. 
For the purpose of calculations it has been assumed that 
and denote maximum and minimum principal strains so 
that 

been 
that 
with 
as a 


q 

the ratio e /e can vary from - 00 to 1 . 

Predictions concerning isoclinic angles have 
shown in Figure 5.1. Since approximate law predicts 
the direction of principal birefringence coincides 
that of principal strain, a plot of isoclinic angle 
function of orientation of principal strain is a 
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straight line at !+5 degree. This is true for all materials 
and any state of strain. However, isoclinic angle predic- 
ted by exact law depends upon both the material properties 
and state of strain. Isoclinic angles for two materials 
subjected to two widely different biaxial states of strain 
have been shown in Figure 5*1* It may be noted that the 
predictions of exact and approximate laws are very close 
to each other for the present material. The predictions 
for the material (Ref. 27) are not as good but they appear 
to be within permissible limits of experimental error. 

Orientation of principal birefringence has been 
compared with the orientation of principal stress for 
different biaxial states of strain in Figure 5.2. In 
this figure all angles were measured with reference to 
the direction of maximum principal strain. As observed 
in Figure 5*1, the predictions of exact and approximate 
laws for principal birefringence direction are quite close 
to each other. It may be noted that the direction of 
principal stress is far away from the directions of prin- 
cipal birefringence or principal strain. A similar conclu- 
sion was drawn for uniaxial state of stress in a study 
by investigators f 18 1 . This conclusion partially jus- 
tifies the development of approximate strain-optic law 
and not the stress-optic law. 
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Predictions of the two laws regarding isochromatic 
fringe order have been sho\>n in Figures 5*3 and 5*4 for 
four different states of strain. The approximate law pre- 
dicts that the fringe order depends only on the difference 
in principal strains and not the orientation of principal 
strain with respect to the fibre direction. Variation in 
fringe order with fibre orientation is predicted by the 
exact law. It is observed that in each case the maximum 
difference between the predictions of the two laws is less 
than 10 percent. In view of the simplification that is 
obtained through the use of approximate law, the error may 
be permissible, particularly for preliminary stress analysis 
problems. 

5.4 ERROR ESTIMATION, CORRECTION AND 
ITS EXPERIMENTAL VERIFICATION 

From the discussion in the thesis so far, it is 
quite clear that the photoelastic response (isochromatic 
fringe order and isoclinic angle) of an orthotropic material 
to a known state of strain can be predicted using exact 
and approximate strain-optic laws. It was observed in the 
previous section that the values of isochromatic fringe 
order and isoclinic angle predicted by two laws do not, 
in general, coincide. The difference arises because, it 
is implicit in the approximate law that the photoelastic 
response is independent of directions of principal strains. 
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Exact theory 


Fibre orientation , a (degree) 


Figure 5.3 Prediction of birefringence from exact and 

approximate theories for principal strain 
ratios of -0.25 and + 0.25 
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It is, therefore, expected that when approximate strain- 
optic law is used to predict magnitude and directions of 
principal strains from the measured values of fringe order 
and isoclinic angle, both the magnitude and directions of 
principal strains will be in error. However, it has been 
observed in the previous section that the maximum error 
in the directions of principal strains is limited to only 
about 7 degree and hence a correction for the same may 
not be essential. It will, therefore, be very desirable 
if a method of estimating error in the magnitudes of prin- 
cipal strains can be evolved. 

It was intuitively felt that the exact strain- 
optic law and a state of strain predicted by the approxi- 
mate strain-optic law may be combined to extract some 
useful information. In particular, it was hypothesised 
that the ratio of exact fringe order to approximate fringe 
order (calculated by substituting approximate strains in 
the exact strain-optic law) is same as the ratio of exact 
and approximate principal strain-difference. Experimental 
verification of the hypothesis along with other details 
and the procedure to be followed in the stress analysis 
problems has been discussed in this section. 

A procedure of determining exact principal 
st rain- differe nee according to the hypothesis, qualita- 
tively stated above, really consists of the following 
steps? 
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It is, therefore, expected that when approximate strain- 
optic law is used to predict magnitude and directions of 
principal strains from the measured values of fringe order 
and isoclinic angle, loth the magnitude and directions of 
principal strains will be in error. However, it has been 
observed in the previous section that the maximum error 
in the directions of principal strains is limited to only 
about 5 degree and hence a correction for the same may 
not be essential. It will, therefore, be very desirable 
if a method of estimating error in the magnitudes of prin- 
cipal strains can be evolved. 

It was intuitively felt that the exact strain- 
optic law and a state of strain predicted by the approxi- 
mate strain-optic law may be combined to extract some 
useful information. In particular, it was hypothesised 
that the ratio of exact fringe order to approximate fringe 
order (calculated by substituting approximate strains in 
the exact strain-optic law) is same as the ratio of exact 
and approximate principal strain-difference. Experimental 
verification of the hypothesis along with other details 
and the procedure to be followed in the stress analysis 
problems has been discussed in this section. 

A procedure of determining exact principal 
strain-difference according to the hypothesis, qualita- 
tively stated above, really consists of the following 
steps s 
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(i) Measure isochromatic fringe order, N, and isoclinic 
angle, 0, at a point in the model for a known load. 


(ii) Using approximate strain-optic law (Eq. 5.10), 

approximate principal strain-difference, ( e - e „) 

5 P q ap 

can he obtained as : 


( e - e ) 
p q ap 


= N f. 


(5.1 ^ 


j>j 


where N has been measured in step (i) and f e 

is the strain-fringe value as obtained in an earlier 

section. 


(iii) According to approximate strain-optic law the direc- 
tion of maj or principal strain coincides with the 
isoclinic so that; 

0 = 0 (5-l4) 

dp 

(iv) The exact strain-optic law (Eq. 4.27) can be written 
in the following form; 


N 


£ L 


{ t ( e p - e 


q f 


£L ) - sin 2 0 


(e - e ) (1 + 

P 1 f eT 


eT 

2 


) ] 


1/2 


+ t (e -e ) 

P 1 


sm 


2 0 ] > 


(5.15) 


It is noted that in Eq. (5.15) the principal strains 
appear only as the difference except in one term where 
e^ is multiplied by the ratio f £ j/i’ e j * 4s an approxima- 
tion this term may also be assumed to be equal to the 
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principal strain— difference. The error introduced by this 
assumption will be small because the ratio f eL /f eT is close 
to 1. Therefore Eq. (5.1 5) can be written ass 

e -e 9 f 2 

IT = -£ {[1 - sin^ 9 (1 + ] 

£L 1 el 

+ [ sin 2 9 ] 2 } 1/2 (5.16) 

eLT 

It may be pointed out that if the ratios and 

f el/ f eL rp are a PP roximated ' to 1» E q« (5.16) will reduce 
to the approximate strain-optic law. By retaining the 
ratios f eL /f eT and f e L /f ^ at all places except one, 

Eq. (5.16) is still quite close to the exact strain-optic 
law. 


>) Calculate approximate fringe order, N , by substituting 

(e -e ) and 9 ^ in Eq. (5«l6). 

P q ap ap 1 

’vi) Obtain the corrected principal strain-difference, 


( e - e ) , from the following relation: 

v p q cr 


( e - e ) 
p q cr 


N 


N 


ap 


( e - e ) 
p q ap 


(5.17) 


Experimental verification of the hypothesis, that the 

( e - e ) is closer to the exact principal strain 
v p q cr 

difference, ( - e n ) . has been accomplished by con- 

ducting experiments on models subjected to uniaxial and 
biaxial states of stress-ydescribed in the following 
paragraphs. 
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Photoelastic behaviour under uniaxial state of 
stress has been studied using tension specimens. Fibre 
orientations with respect to a longitudinal edge of the 
specimen changed and kept at 0, 15, 30, 45, 60, 75 and 90 
degree in order to obtain different states of strain for a 
fixed state of stress. In each case the isochromatic 
fringe order was measured with increasing load. A typical 
curve of fringe order as a function of applied stress has 
been shown in Figure 4.1 and was discussed earlier. The 
isoclinic angle was measured when the effect of initial 
birefringence could be neglected. The principal strains 
and their directions were determined by fixing strain 
gauges on the specimens. The strains can also be calcula- 
ted from the known state of stress and the elastic cons- 
tants of the material. Both the methods gave identical 
results. Measured values of fringe order, isoclinic angle, 
principal strain-difference and the direction of major 
principal strain have been given in Table 5.1 for a fixed 
state of stress (this state of stress was kept the same 
for all fibre orientations). Also given in Table 5.1 are 
the calculated values of ( - e q_) ap j p “ e q^cr and 

errors in the uncorrected and corrected values of principal 
strain-dif f erence . It may be observed that the exact value 
of principal strain- difference is closer to the corrected 
principal strain-difference than the uncorrected one. 


Table 5,1 Results of Photoelastic Analysis for a Uniaxial State of Stress 
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The last two columns show that the maximum error in 
uncorrected principal strain-difference (obtained from 
approximate strain— optic law, Eq. 5*10) is 12 percent where 
as the error in corrected principal strain- difference is 
limited to less than 5 percent. The absolute magnitudes 
of error for the corrected and uncorrected strains have been 
presented in Figure 5*5* This clearly demonstrates that 
the correction indicated by Eq. (5*17) is in good agreement 
with experimental results. 

Photoelastic behaviour under biaxial state of 
stress has been studied using a circular disc under dia- 
metral compression. The observations have been made for 
the cases when the load axis coincides with the fibre direc- 
tion and when the two are at 30 degree to each other. In 
each case measurements concerning fringe order, isoclinic 
angle and state of strain have been made at three different 
points, whose locations are schematically shown in Figure 
5*6. In both the cases the locations of the points with 
respect to the centre of the disc and fibre direction were 
the same. The directions of principal strains and principal 
birefringence (as indicated by isoclinics) have been shown 
in Figures 5.7 and 5*8. Photographs in Figure 5*9 show 
typical isochromatic and isoclinic fringe patterns. The 
actual state of strain was obtained by strain gauge measure- 
ments at the three points. The results of different 
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measurements and calculations have been given in Table 5.2. 
It is cnee again observed that in each case the exact 
principal strain-differenco is closer to the correct- d 
values of principal strain-difference than the uncorrected 
one. It is observed that the percent error for corrected 
principal strain- difference is always less than that for 
uncorrected one. The maximum error in uncorrected principal 
strain difl erence is about + 6 percent whereas for the 
corrected one it is limited to + 3 percent only. This 
further demonstrates that the correction factor indicated 
by Eq. (5.17) is in good agreement with the experimental 
results for biaxial state of stress as well. 

5.5 CONCLUDING MARKS 

An approximate strain-optic law has been deve- 
loped in this chapter which can be easily applied to obtain 
principal strain-difference and the directions of principal 
strains from the photoelastic observations. It has been 
experimentally established that the error in predicting 
the principal strain directions is limited to few degrees 
only whereas the error in the principal strain- difference 
is less than + 15 percent. In view of the simplicity of 
the strain-optic law it is suggested that for preliminary 
stress analysis this error may be permitted. A correction 
factor has been developed which can be applied for a more 
accurate analysis. It has been experimentally established 



Table 5*2 Results of Photoelastic Analysis for a Biaxial State of Stress 
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** Calculated parameters 
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that the correction factor indicated by Eq. (5.17) bri ng S 
?; the corrected value of principal strain-difference within +5 

i- percent of the exact principal strain- difference. 
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Ik 

% 

i, 

l 

lb 

v 




CHAPTER 6 


CONCLUSIONS 

Duo to increasing importance of glass-fibre 
reinforced plastics in various industries, the present 
investigations wore planned to obtain better understand- 
ing of the brittle coating method and the technique of 
transmission photoelasticity and to develop them to a 
state where they can bo easily annlied to the problems of 
stress analysis involving orthotropic materials. Based on 
the studies reported in this thesis the following conclu- 
sions can be drawn. 

1. When a brittle coating analysis is carried out on the 
orthotropic materials, the coating cracks represent 
the directions of the principal strains in the speci- 
men material underneath and not the directions of the 
principal stresses. This conclusion is valid for 
uniaxial as well as biaxial stress fields. Further, 
the threshold strain of the coating is not affected by 
the biaxiality of the stresses in the coating. 

2. Optically superior photoelastic composite materials 
can be obtained by fulfilling the following two condi- 
tions : 

(a) The refractive indices of the constituents of 
the composite should be matched as closely as 

possible. 
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(b) The entrapped air in the composite should be 
minimised. 

3. A poor matching of the refractive indices results into 

the following undesirable effects: 

(a) A parallel beam of light passing through the model 

smeared in the direction perpendicular to 
the fibres. 

(b) Sharp boundaries of an object, when viewed through 
the model 5 appear smeared, again, in the direction 
perpendicular to the fibres. As a consequence 
very closely spaced straight parallel lines will 
not appear distinct. 

(c) White light passing through the model will be split 
in its colours and exhibits a double spectrum 
symmetrical about the centre line. 

The above effects are undesirable for photoelastfc 
applications. The first two cause a reduction in the sharp- 
ness of isoclinic and iso chroma tic fringes. The third one 
will interfere in the analysis if white light is used. 

It is suggested that the matching of the refrac- 
tive indices be judged by the absence of above effects and 
not by visual inspection for transparency. 
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4. An experimental assessment of the exact theories of 
photo-orthotropic-elasticity (including a strain-optic 
law derived in the present work) shows that the photo- 
elastic response of orthotropic materials can he quite 
accurately predicted through these theories, however, 
no theory permits, evaluation of the individual strain 
or stiess components or even the difference between 
their principal values from birefringence measured in 
a polari scope. The directions of principal stresses or 
strains also can not be established. 

5. An approximate strain-optic law has been derived from an 
exact strain-optic law which in turn was derived from 
the well established theory of photo-elasticity of 
crystals. This law can be easily applied to obtain 
principal strains from the photoelastic observations. 

It has been experimentally established that the error 
in predicting the principal strain directions is limited 
to few degree only, and the error in principal strain- 
diffcreneo is less than + 15 per cent. 

6. A method has been developed by which a correction can 
be made in the approximate value of principal strain- 
difference. It has been experimentally established 
that the corrected value of the principal strain 
difference is within + 5 per cent of tlle exact P* 11101 ” 


pal strain-difference*. 
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APPENDIX A 


THAN SFOEM ATIO N RELATIONS FOR STRESSES 
AMD IMATEKEAL PROPERTIES ’ 


STRAINS 


For an urtholropie iamina 3n the 1-2 plane (or L-T plane 
where I. and T are the longitudinal and transverse symmetry 
ax ca of blit' material) as shown in Figure A.1, the strain- 
stress relation, in berms ox engineering constants, for a 
plane stress state may bo written as [ 35 ] : 
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If an a;xln nonnal to 1 -I! plane is designated as 3, the 
Ecu (Ini) neodn to bo supplemented by the following equa< 

tions t 
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The si raiu-sb rose relation as given by Eq. (A.1) is 
referred to the principal material directions. However, 
the coordinate axes chosen for problems often do not 
coincide with the principal directions of ortho tropy. 
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necessary to transform the strain-stress 


lb i s , therefore, 
rrnntloMi; fro,,' one coordinate system bo another. 

l )0 b IJimv in? n,no Uior X - Y coordinate system 
makinr, nn the 1 or L - axis, the transfoma- 

tioii:; for Un 1 s trasses and strains can be written as; 
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ls a transformation matrix given as: 
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, . /, n in X - I coordinate 

The strain-stress relation, C • 
system iss 
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And the transformed engineering constants in terms of 
principal material -properties can be given ass 


(A.6) 
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appendix b 


failube theory 

In order to predict stresses in the specimen from 
observations of the crack patterns In the coating, some 
coating failure law must be utilised. Three laws of failure 
maximum strain theory, maximum stress theory and the Mohr 
theory - have in the past been considered to represent the 
behaviour of coating. The maximum- strain and the m aximu m 
stress theories are not verified exp eriment ally. The Mohr 
theory of failure is frequently employed to predict failure 
stresses in materials for which the ultimate compressive 
strength ( a uc ) is much greater than the ultimate tensile 
strength ( ) . Since it is well known that the compre- 

ssive strength of a brittle coating grfeatly exceeds its 
tensile strength, the Mohr theory of failure would be appro- 
priate for predicting its failure. 

In this application of the Mohr theory of failure, 
three important cases arise; 

c c 

(i) For biaxial tension stress state ( CT p > a q > 0)» 
failure of the coating occurs when 



For this case the predictions of Mohr theory 
coincide with those of maximum stress theory of 


failure . 


(ii) 
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For tension - compression state ( c c 
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> 0 > a 


the failure relationship becomes 
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(iii) 


For biaxial compression state (0 > a p C > a^) the 
failure equation will be 
„ c • 


-> 1 
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uc 


( B a 3 ) 


In this formulation it is usual to assume that 
_ c . 

= K (BA) 


uc 
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ut 


In order to obtain the coating stresses from these 
failure equations, it will be necessary to evaluate the 
tensile and compressive strengths of the coating, intimate 
tensile strength of the coating is obtained by applying 
coating on a cantilever beam of "Hie specimen material. The 
stresses in the beam under an end deflection are : 
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The minimum strain at which a crack appears at 
any point is given by = e where, £ is the 

threshold strain of the coating. Now introducing these 
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conditions in Eqs. (2.7) and solving for coating stresses. 


the ultimate tensile strength of coating is obtained as 
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The stress a ' is given as 
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It may be noted that uniaxial stress in the specimen pro- 
duces biaxial stresses in the coating due to mismatch bet- 
ween Poisson’ s ratios of the coating and the specimen. It 
is established practice to calculate from Eq. (E.6) 

because the value of o q ° given by (B.-7) is comparatively 
very small. Further, it has been established in Chapter - 2 
that the threshold strain, e^, is not affected by the biaxia. 
lity of stresses in the coating. 



APPENDIX C 


EVALUATION OF STRAIN - OPTIC PROPERTY (P 22 - P^) 

As discussed in Chapter 5, the 2-3 plane of 
transversely isotropic composite behaves macro scopically 
as an isotropic plane to in-plane loads so that a theory 
applicable to mi croscopi calls'- heterogeneous but macrosco- 
pically isotropic material can be applied to obtain elastic 
and photoelastic properties of such a composite, Richard 
and Young [39] developed a composite which behaves in an 
isotropic, homogeneous manner on a macroscopic scale. The 
constituents of this composite are solid-glass microspheres 
embedded in a polyester matrix. The composites with varia- 
ble mechanical and optical properties were obtained by 
changing the volume fractions of the constituents. For 
predicting the photoelastic properties of such composites, 
a theory has also been developed by them. According to 
this theory, the stress-fringe value f of the isotropic 
composite can be given as; 
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where V g is glass content; f ag and f 0m are the stress- 
fringe values for the glass and matrix respectively, and 
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and X m are respectively the ratios of principal stress 
difference in glass and matrix to that in 'the composite. 
The values of X and X^ can he given as: 
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In these expressions the average coordinate stresses 0 

a 0 ™ and a , m can b® evaluated by the following 
yjjj yxu 

relations : 
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Evaluation of average coordinate stresses given by Eq. (C.3) 
will require the values of elastic constants of the consti- 
tuents together with the transverse modulus Ej and Poisson's 
ratio, v 2 3 ? of the composite. Due to practical difficul- 
ties associated with the direct determination of it 

was decided to calculate it theoretically following 
Bichard [4-0] . Bichard has shown that the mechanical pro- 
perties of glass-micro sphere-filled composite can be 
evaluated with a reasonable degree of accuracy using a 
self-consistent approach proposed by Hill [42, 431. The 
properties of constituents for uni-directional model 
material developed in the present work are given in Table 
C.1, Knowing the constituents properties, volume fraction 
of glass fibres and the elastic properties Erp and v £3’ 
the stress-fringe value, f CTG , f° r a 2-3 plane-composite 
(light passes normal to 2-3 plane) can be evaluated from 
Eq. (C.1). The strain-optic property, (P 22 ~ P 23^ } can " be 
calculated from the following expression: 

Em 


(P 22 " p 23 } 
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Substitution of the Eq. (C.1) into (C.5) will yield: 
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1 + v 23 u ° “eg 
The calculated properties, v 2 2 and (1*22 “ ^23^ ^ or 
rectional material (v„ = 0.165) were found to be 0.34 and 


0.1016 x 10° fr/m respectively. 
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Table C,1 Properties of the Constituents of 
Unidirectional Model Material 


Material 

[^Stress-fringe J 
[ Value 1 

1 10 3 N/f r~m ! 

Modulus of 
Elasticity 

10 6 N/m 2 

Poissons 1 
Ratio 

Fibres 

(E-glass) 


238.2 

72408.0 

0.20 

Matrix 
('Polyester 
+Dibutylphthala to ) 

32.5 

4172.0 

O .36 

*with X = 

5893 A 
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